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INTRODUCTION 


This  final  report  Is  part  of  a continuing  effort  to  study  and  under- 
stand free  radical  processes  of  importance  in  combustion  and  other  chemical 
systems . 


We  include  two  preprints  of  papers,  Chapters  1 and  2,  submitted  to 
the  Journal  of  the  American  Chemical  Society.  These  papers  describe  our 
recent  work  on  combination  reactions  and  heats  of  formation  of  stabilized 
radicals . 


Chapter  3 is  a preprint  of  a paper  submitted  to  the  International 
Journal  of  Chemical  Kinetics.  This  paper  is  complementary  to  our  initial 
work  on  multiphoton  dissociation,  which  was  published  (J.  Amer.  Chem.  Soc . , 
99,  8063  (1977)  ) and  is  included  as  Chapter  4) . 


Chapter  5 represents  work  of  both  practical  importance  for  stratos- 
pheric modeling,  as  well  as  theoretical  Importance  with  respect  to  under- 
standing the  details  of  radical-radical  Interactions.  This  paper  recently 
appeared  in  the  International  Journal  of  Chemical  Kinetics,  10 , 489  (1978)  . 

Chapter  6 has  been  submitted  to  the  Journal  of  Chemical  Physics.  This 
paper  illustrates  our  radical  combination  model  for  the  prototypical  methyl 
radical  combination. 
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I INTRODUCTION 


Allyl  radical  is  the  prototype  of  a resonance  stabilized  radical  - 

thermochemistry  and  reactivity  in  the  gas  phase  has  been  the  subject  et 
numerous  investigations.*"^  There  has  been  considerable  controversy**  -*»-- 
the  correct  value  of  the  heat  of  formation  and  the  allyl  resonance  eeer* . 

(ARE)  Quoted  values  for  ARE  range  from  10  to  25  kcal/mol,  but  recent 
experimental  values  for  ARE  now  seem  to  fall  around  11+2  kcal/nol  . Tb« 
correct  value  for  ARE  is  certainly  fundamental  for  a thorough  under standii^ 
of  chemical  bonding  and  ground  state  properties  of  conjugated  radicals. 
Theoretical  calculations  on  open-shell  species,  such  as  allyl  radical,  a»« 
hampered  somewhat  by  intrinsic  difficulties,®  but  despite  this  problnn.  «be 
Generalized  Valence-Bond  (GVB)  concept  appears  to  be  quite  successful  .* 

Furthermore,  the  need  for  the  accurate  determination  of  rndical-radica . 
recombination  rate  constants  has  been  clearly  pointed  out.^  Absolut*  rst* 
constants  for  many  disproportionation  reactions  and  radical-molecule  rmactiam* 
are  critically  dependent  on  the  rates  for  combination  of  the  radicals,  sisce 
many  rate  constants  have  been  measured  relative  to  the  radical  combinatios 
rate  constants.  Thermochemical  parameters  of  free  radicals  can  also  be 
obtained  from  the  Arrhenius  parameters  of  free  radical  combination  in  cases 
where  the  reverse  reaction  (dissociation)  has  been  studied.  Radical  recom- 
bination rate  studies  in  the  case  of  resonance  stabilized  radicals  are 
sparse  and  the  rates  for  allyl  and  2-methallyl  recombination  have  only  bees 
measured  at  ambient  temperature  by  flash  photolysis.®*®  It  seemed  approprist* 
to  extend  the  rate  measurements  in  order  to  get  a reliable  set  of  Arrbeaiue 
parameters  for  allyl  recombination  at  higher  temperatures  where  this 
reaction  is  the  prototype  for  important  chain  terminations. 
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In  this  paper,  we  report:  j 

(A)  The  equilibrium  and  kinetics  of  allyl  radical 
recombination  at  844  < T/K  < 1061,  and 

(B)  The  recombination  kinetics  of  allyl  radical  at 
T = 625  K. 

The  method  used  is  Very  Low-Pressure  Pyrolysis  (VLPP)  employing  a newly  j 

designed  molecular  beam  sampling  apparatus  thereby  eliminating  complicated  i 

secondary  reactions  of  radicals  on  the  walls  of  the  mass  spectrometry  ^ 

chamber 


H KXrKKlMKNTAl, 


KlKur«  1 dt*plct8  the  Koneral  uxpurlmentn I duHlKit  of  the  moleoiilnr 

heum-HampllitK  VIJ^P  itppurntuit.  A vi« I vo-capl  1 fury  tub«*  a rruiiKoniont  deliver!* 

a contix>llod,  steady  flow  of  reactant  fo  the  reactor  (typically, 

10**  - S’*  10***  molocnleN.  aec)  . In  cases  where  the  reactant  kun  bad  an 

vmsultahly  low  vapor  preasure  at  room  temperature,  the  whole  tins  Inlet 

system  was  placed  In  a hot  hox . The  flow  then  enters  the  low-pressure 

reactor,  and  subsequently  effuses  from  t!>e  reactor  exit  aperture.  The 

product  nas  flux  Is  collimated  (<  I cm  In  diameter)  on  passing;  tl\ri>UKh  a 

differential  pumpli^  cliamber,  modulated  with  a rotatluK  chopper  wheel, 

and  finally  lont/.ed  and  detecteil  by  a quadrupole  mass  spectrometer 

(Ktnnl){an  -too).  A lock-tn  amplifier  (Princeton  Applied  Itesearch  Model 

1284)  scv>arates  the  modulated  slKual  from  tl\e  unmodvilated  backKCound 

slKUal.  Kluure  2 dtsytlays  the  two-av>erture  reactor  deslRn.  Ueactor 

parameters  are  Riven  by  (It  ■ lai'RO  aperture,  K ■ small  aperture): 

V .134  I,  M M 1082  ^ (T  M)‘  » S-*  , k ^’(11)  2.S.S71  " (T  V)*  “ s“*  , 

e 

k (S)  ■ .2088  ' (TV)*  * s"*  . The  all-qunt'f^  reactor  (t.e.,  Knudsen 
e 

cell)  encased  In  a nickel  block  was  heated  with  an  I'lectrlcal 
clam-shell  heater,  and  the  temperature  was  measured  by  a chromt'l- 
alumel  thermocouple.  The  experiments  were  carrleil  t>ut  by  monllortitK 
mass  sped  romet  rl  c intensities  as  a function  of  th»'  t empera  t tire , 
the  flow  rate  of  the  Ras  Into  the  reactor  (K^^  molecules  s“*)  and  th«*  resi- 
dence time  of  the  molecules  within  the  reactor  (small  or  laiye  aperture).** 
niallyloxalate  (03ll,^0t'0t'CKX'.,ll^ ) was  purchased  from  PfaHv.  and  Hauer,  Inc., 
and  Was  purified  simply  by  pumpliiR  off  all  lower  bollliiR  fractions  (mainly 
l'll3‘=Cll-l'H3-('llO)  at  10"*  torr.  IK'-VS-analysls  revealed  no  further  Imvnirt- 
tles  present.  3 ,3 '-a/o-l-propone  (I'jIl^N  ) was  prev>ared  accordtiiR  to 
literature  procedures.*'*  (The  samples  contained  H|,0  and  dlethylether 
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which  did  not  interfere  with  our  measurements.  The  concentration  of 
3 ,3 ^-a^o-l-propene  was  measured  by  monitoring  m/e  = 28  in  the  absence  of 
the  molecular  ion  of  at  m/e  = 110.) 

The  reactions  were  all  studied  in  a reaction  vessel  which  was  "cleaned" 
after  every  set  of  two  or  three  experiments  (flow  rate  studies)  by  passing 
air  through  the  vessel  at  temperatures  around  1100  K.  This  procedure  served 
to  remove  the  soot  which  was  formed  at  high  temperatures  (1100  K)  and  could 
be  seen  on  the  walls  of  the  reaction  vessel . Such  a procedure  was  necessary 
because  the  carbon  coating  greatly  enhanced  the  interference  by  heterogeneous 
reactions . 
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III  RESULTS  AND  DISCUSSION 

A.  Equilibrium  and  Kinetics  (844  ^ T/K  ^ 1061) 

In  the  range  844  ^ T/K  < 1061,  the  following  reaction  system  was  studied 
using  diallyloxalate , CgHj^O^ (DAO) , as  a convenient  source  of  allyl  radicals:*’ 


C3H5OCOCOOC3HS (DAO) 

■*‘2.  > 2C3H3-  + 2CO3 

(1) 

aCaHj* 

^*P^CeH,„(BA);  Kr,d  = V^d 

(2) 

DAO 

k 

C3H5OCOCOOC3H5 

(3) 

k 

CaHs- 

e 

(4) 

kBA 

^6**10 

G 1 

(3) 

CO2  served  as  an  internal  standard  or  monitor  for  allyl  radical,  because  for 
every  allyl  radical,  a molecule  of  CO^  was  formed.  This  was  verified  exper- 
imentally by  "titrating"  CgHg*  with  DI : 

CjHs*  + DI  2 ►CjHjD  + I’  (6) 


For  [Dl]gg  ^ 5 * 10"®  M,  it  was  verified  that  [CaH^Dlgg  = [cO^lgg  through 
independent  calibrations  with  propylene  and  CO, • Furthermore,  this  is  an 
experimental  demonstration  that  the  formation  of  1 ,5-hexadiene  proceeds  through 
recombination  of  Intermediate  allyl  radicals  and  not  through  unlmolecular 
elimination  from  DAO.  It  was  found  that  the  disproportionation  reaction  (7) 
is  not  Important  under  our  experimental  conditions: 


aCjHs-  + C3H,, 


(7) 
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a conclusion  further  confirmed  by  the  absence  of  sizable  amounts  of  allene^ 
and  by  the  results  of  the  Braumnn  plot  (Vide  Infra)  . 

The  steady-state  kinetic  expression  for  the  reaction  system  (1)  to  (5) 
In  n low-pressure  stirred  flow  reactor  Is  (see  Appendix  for  details  of  the 
derivation)  : 


k + k 
d e 


BA 


BA 

= k * 


2R" 

BA 


- - k - 

- 2R‘>  )“  e 
BA 


(8) 


BA  r f V ^ « 

where  k Is  the  escape  rate  constant  (In  s"  ) for  1 ,5-hexadiene  (BA), 
e 

R"  Is  the  specific  flow  rate  (In  moles  s”*  I”*)  of  BA  out  of  the  reactor 
BA 

(as  monitored  by  m/e  = 67  or  54)*^  and  R®  Is  the  corresponding  flow  of 

COj 

COj,  (as  monitored  by  m/e  = 44)  Two  experiments,  each  with  different 
BA  , 

k (or  residence  time)  were  performed  at  a given  temperature,  varying 

f 

K typically  from  10‘*  - 4 10*®  molec  s“*  , so  that  k and  k could  be 

determined  separately,  but  this  was  possible  only  over  a limited  temperature 
range.  Table  I displays  the  results  obtained  over  the  temperature  range 
(844-1061  K)  and  Figure  3 shows  the  corresponding  van't  Hoff  plot  for  K 

r.d 

The  results  can  be  summarized  as  follows: 


k 

r 


ln(K  ^/M"‘) 
r.d 


(1.9  + .8)  * 10®  M~‘  s"‘ 

-33.5+  10.71  54540 

R ^ RT  ’ 


<T>  = 880  K 

<T>  = 950  K 


(9)»» 


Thus  this  method  provides  a direct  measurement  of  equilibrium  (£)  with 

AS®  = -33.5  eu,  AE®  = -54.5  kcal/mol  and  AH®  = -56,4  kcal/mol.  The  variable 

and  sometimes  large  error  limits  of  K displayed  In  Figure  3 are  the  result 

r ,d 

of  the  algebraic  separation  of  k and  k from  two  Independent  experimental  sets 

r d 

DA 

of  k (large  and  small  apertures)  according  to  (8).*®  The  heat  capacity  data 
e 

for  1 ,5-hexadlene*^  and  allyl  radical*®  are  such  that  AS®  and  AH®  show  no  tem- 
perature dependence  from  300  to  930  K.  With  AH^(BA)  = 20.2  kcal/mol  (Table  II), 
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and  = —56.4  ±.1.0  kcal/mol  (.figure  3),  one  obtains  ^f(allyl,g)  - 38.3  ±1.5 
kcal/mol  at  T = 300  K.*®  The  second  law  e.itropy,  AS®  = — 33.5Je.u.  ( (9),  .stan  - 
dard irate:  1 atm)  compares  favorably  with  an  a priori  estimate  of  —34.80  ±1.0 
e.u..*’**“  T.he  present  result  Indicates  an  allyl  resonance  energy  (ARE)  of 
12.5  kcal/mol  on  the  basis  of  a value  of  98.0  kcal/mol  for  BDE  of  propane. 

This  is  in  good  agreement  with  the  reported  literature  values  from  shock-tube 
decomposition  studies  of  olefins^ ^ and  in  excellent  agreement  with  a value 
determined  by  photoionization  mass  spectrometry^  of  38 .4  ± 1 .7  kcal/mol  for 
^j(allyl,g).  The  agreement  with  a value  of  39.10  +1.0  kcal/mol  for  AH^ 
derived  from  the  measurement  of  the  equilibrium  (10) 


/\  + 1- 


A. 


is  also  noteworthy. 


It  was  pointed  out  in  the  Introduction  that  the  earlier  study*  of 
the  present  system  (l)-(5)  yielded  values  for  and  kj.,  somewhat  too  high, 

thus  providing  a value  of  AH?  (allyl, g)  = 41.2  kcal/mol  and  ARE  =9.6 

X y 3 0 O 

kcal/mol,  together  with  k^  = 6.2  + 10®  M“*  s~*  at  913  K Lj  1033  K.  The  reason 
for  this  discrepancy  was  found  to  lie  in  a wall-catalyzed  recombination 
reaction  in  the  mass  spectrometry  chamber.*®  The  present  experimental  con- 
figuration of  the  apparatus,  however,  rules  out  complicating  reactions  of  this 
kind.  The  apparent  excellent  agreement  with  GVB  calculations®  should  be 
taken  "cum  grano  salts,"  because  the  "stabilization  energy"  is  referred 


to  the  total  energy  of  the  hypothetical  canonical  structure 


, which 


is  not  an  observable  species  and  therefore  experimentally  not  accessible. 

The  measured  value  for  k^.  (9)  is  discussed  together  with  values  obtained 
at  T = 625  K in  section  B . 


A 


i 


Table  11 


THERMOCHEMICAL  QUANTITIES 


1 ,S-hexudlene 


allyl  radical 


kcal/mol 


Table  III 


RATE  CCWSTANTS  FOR  ALLYL  RADICAL 
RECOMBINATION  AT  T = 625  K 


n 

qO 

^3  00 

C “ 

P ■ 3 0 U 

a 

Co 

8n.40 

62  .10 

28  .80 

14  .60 

50.10 

24.95 

B.  Keeomblnat Ion  Kinetlc.s 


4 


1 


The  recombination  kinetics  of  allyl  raitlcul  at  T - 625  K were  stuillea 
usliiK  3 ,3  *-aa!o-I-propene  (A^)  as  a source  of  allyl  radicals.  The  following 
reaction  scheme  was  u.ssumtsi; 


k.. 


kjj^  j*  ^ * 
> 


111) 

1 13) 

113) 

111) 


The  steady-state  kinetic  expression  Ithe  Urauman  equation)  for  the  reaction 


scheme  (ll)-114)  In  a sttrred-flow  reactor  Is: 


i t 


.\l. 

K“ 

lU 


(k 


I + 


k )V 
w 


l/s 


UR) 


3kr 


Vt  |,<o  IJ 
lU 


where  V'\  Is  the  flow  (molecules  s"*)  of  3,3  -a/o-l -propene  Into  the  reactor, 
AZ 

K**  Is  the  flow  of  1 ,R-hexadlene  cut  of  the  reactor  (as  monltor»>d  at  m e ■ 

C 11  • 

\;7  or  .M)  , k^,  ^ * Is  the  escape  rate  constant  (s“‘)  of  allyl  radical  and 
k ls“*)  Is  the  rate  constant  for  a first-order  loss  process  of  the  allyl 

’V 

radicals  (nuist  likely  on  the  walls)  . I'he  temperature  of  1>2R  K was  chosen 

such  that  A/,  was  decomposed  almost  quantitatively,  so  that  the  mass  spec- 

trometrlc  Intensity  of  me  = 38  (corrected  for  the  contribution  of 

I ,5-hexadlene)  reflected  the  amount  of  allyl  radicals  found  In  the  reactor 

(K^  ) after  assessment  of  the  calibration  constant,  >•1  (ustn*:  air),  which 
Az  ^ 

relates  the  MS  Intensity  (m  e * 38)  to  the  rate  of  production  (K  » K*)  ) or 
(steady  state)  concentration  of 
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A!\  alternative  procedure  for  determination  of  c.>n!ilatt>d  of  "titration" 

A.’. 

of  with  01  it*),  while  monitoring  the  MS  Intensity  of  m/e  » -13  ICjH^O) 

under  the  condition  -*  ®,  so  that  the  amount  C,H^I)  formed  was  repre- 

sentative of  the  amount  of  allyl  radical  orlKlnally  present.  However, 
the  concentration  of  C,Hj,‘  had  to  be  kept  low  I [CjHj  * ] ^ 3.7  >t  10** 

molecules  1“*,  correspondlnjt  to  "X  5 >i  10**  molecules  s“‘  at  T ■ 625  K)  , 

and  Loi]  had  to  be  quite  high  In  view  of  the  smaller  rate  constant  for 
ss 

D-transfer  Ik^),  compared  to  the  recombination  rate  constant  Ik^,)  ; 

kj. /kj  er  6>1  at  T 625  K.*“  The  disappearance  of  HA  In  favor  of  the  Increase 

In  Cj  l^l)  could  readily  be  observed  at  m/e  > 62,  5-1  (for  HA),  and  43  (for 

CjHjO)  . Both  methods  of  assessing  gave  essentially  the  same  result  to 

within  15%,  but  the  first  procedure  was  used  for  all  reported  e.xperlment  s . 

A plot  of  /'V**  Versus  1 /T**  *''*  gave  straight  lines  acconllng  to  (15)  for 
A*  BA  ^ ^BA  e s •• 

two  different  values  of  k^^  * * * through  the  common  Intercept  1.0  (within 

experimental  erx'or,  see  Klgure  4)  and  Indicated  that  allyl  radical  was 

unable  to  abstract  hydrogen  In  a disproportionation  reaction.  In  line  with 

the  results  of  gas-phase  pyivlysls  experiments  In  static  high-pressure 

reaction  vessels*'’  and  the  resxtlt  s of  Section  A (eqn,  7).  However,  It  appears 

that  allyl  radical  does  undergo  wall  reactions,  because  the  ratio  of  the 

slopes  for  the  large  and  snuxll  aperture  la  Klgure  4 Is  considerably  less 
OH*  0*1* 

than  the  ratio  k * * lB)/k^*  ' (S)  Indicating  a non-negllglble  term  k„.(rirst 
e e 

order  loss  process  of  allyl  radical  according  to  (13),  presunxably  taklitg  place 
on  the  walls**).  The  results  of  the  pyrolysis  of  3 ,3  ^-a/.o-l -p»x*pene  are 
displayed  In  Table  11  from  which  It  can  be  seen  that  the  average  value  for 
kj.  • (6.50  i.  I .0)  '*  10*'  M"*  8”*  corresptntds  closely  to  the  one  derived  at 
room  temperature  using  kinetic  flush  sv»ect roscopy*  (k,.  ■ (8.50  ± 3.0)  >t  10** 

M-*  s-*)  . 
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IV  RKKM  CALCU1ATI(»<S  AND  DISCUSSION 


a 


The  rate  constants  for  recombination  of  allyl  radicals  were  calculated 
using  RRKM  theory^*  for  the  unimolecular  decomposition  of  1 ,5-hexadiene 
making  use  of  the  thermodynamic  relation: 


" *^r,d  * '‘d 


(16) 


Two  different  transition-state  models  were  used:  (1)  the  vibrational  and 
(2)  the  rotational  model.*® 


In  the  absence  of  a complete  assignment  of  the  vibrational  spectrum  of 
BA,  the  vibrational  frequencies  of  BA  were  deduced  from  those  of  propylene. 
The  moments  of  Inertia  for  BA  were  computed  using  standard  bond  lengths  and 
bond  angles.  The  vibrational  frequencies  were  then  adjusted  so  that  the 
computed  entropy  matched  the  value  from  group  additivity  data.^*  Table  IV 
displays  the  results  of  the  frequency  assignments  for  BA. 

1 . The  Vibrational  Model 

In  this  model  for  the  transition  state  of  the  bond-breaking  process, 
the  central  -C-C~  bond  was  extended  2.5  times*®  to  3.85  A.  The  torsion 
around  the  central  -C-C—  bond  becomes  a free  rotation  of  both  allyl  radical 
fragments  in  which  resonance  stiffening*®  of  the  former  two  C-C=C  (torsions) 
in  BA  took  place, and  the  external  rotation  around  the  central  -C-C—  bond  was 
made  active  to  allow  these  two  rotational  modes  to  share  in  the  random 
distribution  of  molecular  energy.  Four  C-C  rocking  modes  of  the  C3H5 -units 
in  BA  had  to  be  replaced  by  unusually  low  frequency  bending  modes  (42  cm-*) 
in  order  to  yield  Arrhenius  A-factors  for  decomposition  of  BA  which  would 
correspond  to  the  measured  rate  of  recombination®  of  (8.50  ±.  3.0)  * 10®  M"‘ 
s-‘  at  T = 300  K through  (16),  see  Table  IV.  The  activation  energy  for 
the  bond-breaking  process  (E^j)  was  set  equal  to  -AH®  for  equilibrium  (2) 


Tabic  IV 


MOLECULAR  PARAMETERS  FOR  RRKM-CALCUIATIONS : 
1 ,5-HEXADlENE  (BA) 

(Standard  State:  1 Atm) 


Molecule 

Complex:  "vibration" 

Complex:  "Rotation" 

3000  (6) 

3010  (2) 

3010  (2) 

2960  (2) 

3000  (8) 

3000  (8) 

2850  (2) 

1350  (6) 

1350  (6) 

1650  (2) 

1300  (6) 

1300  (6) 

Frequencies 

1470  (2) 

1100  (6) 

1100  (6) 

1420  (2) 

950  (2) 

950  (2) 

and 

1400  (2) 

400  (4) 

400  (4) 

1300  (2) 

350  (2) 

350  (2) 

degeneracies 

1220  (2) 

42  (4) 

1170  (2) 

920  (3) 

900  (2) 

800  (2) 

350  (2) 

320  (6) 

90  (1) 

70  (1) 

50  (1) 

1 .54 

3.85  (P"^  = 2.5)“ 

3.85  (?"■  = 2.5)“ 

I I (gr  cm^)’ 
ABC 

1.78  * 10’ 

8.60  * 10’ 

8.60  * 10’ 

10‘®  1,  /gr  cn» 

— 

21 .20 

21 .20 

10®®I,I,/(gr  cn»)® 

— 

— 

(91 .44)’ 

♦ 

I /I 

— 

4.14 

4.14 

Ejoo  /kcal/mole 

— 

55.70 

56.70 

— 

56.40 

56.40 

®300  / 

88  .83 

101.04 

102  .91 : T1  = 90  .0% 

— 

15.90 

16.31  : H = 90  .0% 

log(A,  /m"‘  8-‘) 

r,3oo  ' 

T “ * 

10.40*’ 

h 

10.81;  '^^90,0% 

9 ^ , where  r'^  ia  the  distance  of  both  allyl  radical  fragments  In  the  activated  complex 

and  r”  la  the  central  C-C  bond  length  In  the  molecule. 

'’2,303R  log  : = - 33.50  + 8.35  with  AS®  = - 33.50  e.u.  (exp.) 

' ' 18 

i 

• i 

J 


a^sumlnt;  £ero  activation  energy  for  recombination  (Kj.)  at  0 The  re^iults 

of  the  KRKM  calculations  are  dlsplayeil  In  Figure  5 anil  are  listed  In  Table  V, 
together  with  the  hlgh-pressure  Arrhenius  parameters.  It  Is  obvious  from 
Figure  5 that  the  vibrational  model  successfully  fits  the  experimental  data 
In  this  case  of  recombination  of  two  resonance  stablll/.ed  radicals.  Fast 
e.xperience  with  the  recombination  of  simple  alkyl  radicals  such  as  tert-butyl, 

Isopropyl,  and  ethyl  has  Indicated  that  such  a fixed  transition  state  model 
was  not  suited  to  represent  the  loose  transition  state  for  the  bond  fission 
of  a simple  alkane.  The  model  predicts  a transition  state  with  too  large  a 
heat  capacity  at  high  temperatures  due  to  the  loosening  of  the  four  rocking 
modes  In  the  transition  state,  thus  predicting  Increasing  Arrhenius  parameters 
and  Increasing  values  for  kj,  with  Increasing  temperatures.  The  experiments 
Indicate  otherwise,  however.  In  the  present  case  of  resonunce-stabl  llz.ed 
radicals,  the  above  Increase  In  heat  capacity  Is  essentially  balanced  by  the 
resonance  stiffening  of  Internal  rotations,  thereby  reducing  the  heat  capacity 
of  the  transition  state.  The  net  result  Is  a small  decrease  In  heat  capacity, 
and  a concomitant  small  decrease  In  Kj  and  A^  or  a slight  negative  temperature 

CO 

dependence  (see  Table  V).  Although  k^  Is  seen  to  increase  by  a factor  of  4 
from  300  K to  1000  K,  the  vibrational  model  fits  the  experimental  data  Ikj.) 
for  allyl  radical  recombination  quite  satisfactorily  In  contrast  to  normal 
alkyl  radicals,  where  the  high-pressure  e.xperimenta  1 A-factors  for  recom- 
bination or  bond  scission  have  a pronounced  negative  temperature  dependence. 

2 . The  Rotational  Model 

This  Is  essentially  a Gorin  model,  which  represents  the  four  low-frequency 
bending  vibrations  of  the  two  allyl  fragments  as  hindered  rotations.  The 
Internal  modes  of  the  transition  state  are  simply  the  vibrations  and  rotations 
of  the  Independent  fragments.  This  Is  an  alternative  picture  for  the  transition 

I 
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state  typical  of  simple  bond  scission  reactions.  The  hindrance  (Tl/%) 
is  accomplished  In  the  calculation  by  decreasing  the  effective  moment 
of  inertia  of  each  of  the  two  two-dimensional  ally!  rotors: 


I 

L 


I 


eff 


I 100  - T\  \V^ 

\ 100  j 


(17) 


where  1^  and  are  the  two  one-dlmenslonal  component  moments  of  Inertia  of  an 
allyl  fragment  excluding  the  component  of  the  moment  of  Inertia  around  the 
axis  parallel  to  the  bond  being  broken.  The  hindrance  effectively  decreases 
the  number  of  available  rotational  states  by  confining  the  rotational  motion 
of  each  of  the  fragments  • As  the  centrifugal  barrier  moves  to  smaller 
r'*’  values  with  Increasing  temperature,  ^ Is  expected  to  Increase  with  tem- 
perature, thereby^®  reducing  k”.  The  relation  of  T]  to  the  Arrhenius  A-f actor 
for  decomposition  is  simply: 


Ah/A 


100  - T] 
100 


(18) 


where  Ajj  is  the  A-factor  derived  from  a transition  state  model  with  Tf]  0 . 

In  the  case  'jf  the  rotational  model , the  same  molecular  parameters  were  chosen 

as  for  (1),  Ejj  = — Ah° , except  that  the  four  low-frequency  bending  modes  at  4^ 

cni“^  were  replaced  by  two  two-dimensional  rotations  of  the  allyl  fragments 

FdOO  - Tl)!  V'2 

(1  = 8.36  * 10®®  * I (gr  cm®)®  ) . It  was  found  that  a good 

eff  L 100  J 

fit  to  the  experimental  data  could  be  obtained  by  choosing  the  following 

values  for  T];  T = 300  K,  11  = 90%;  T = 625  K,  1]  = 91.5%;  T = 900  K,  11  = 93%. 

Table  V shows  the  results  of  the  calculation  and  Figure  5 displays  the  corres- 

00 

ponding  plot  together  with  the  experimental  rate  constants  k^.  and  k^..  Not  too 
much  Importance  should  be  attached  to  the  absolute  value  of  11,  because  this 

number  Is  dependent  on  the  details  of  the  transition  state  model  with  respect 

CO  Its  structure  and  the  vibrational  frequencies.  Table  V describes. 
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thert'foro,  only  a "reasonable”  choice  of  molecular  parameters  for  transition 
stales  (1)  or  (2).  Furthermore,  llie  exact  relationship  betiveen  H anil  the 
jceometrlcal  parameters  of  the  fragments  (rotational  "freedom"  or  tightness) 

In  the  transition  state  are  virtually  unknown.  It  should  be  noted  that  In 
this  case,  the  dependence  of  H on  tempt-rature  Is  exceptionally  weak  In 
comparison  with  two  recent  examples.^'*  For  methyl  radical  recombination, 

H Increases  from  03', r to  83%  over  the  temperature  range  3l)()-l-100  K and  In  the 
case  of  the  recombination  of  IIO^  ' with  N0._j  • , >1  varies  from  92%  to  98%  over 
the  temperature  range  217-300  K. 

An  Intuitive  choice  for  the  Interfragment  distance  In  the  critical  con- 
figuration of  a normal  alkane  undergoing  bt>nd  fission  to  two  alkyl  radicals 
would  be  the  top  of  the  centrifugal  barrier  .* ’•**’  Kecent  calculations  on 
ethane  dissociation  suggest,  however,  that  the  application  t>f  the  criterion  of 
minimum  density  of  states^*  Is  in«>re  appropriate  In  locating  the  critical 
configuration  with  respect  to  the  reaction  ct'orillnate . if  no  other 
structural  changes  In  the  molecule  occur  during  the  bi>nil  split,  the 
location  of  the  minimum  density  of  slates  will  coincide  with  the  location 
of  the  top  of  the  centrifugal  barrier.  As  the  V.-i'.  bond  ts  breaking,  however, 
lowering  of  the  rocking  modes  Is  expected  to  occur.  ixnisenlng  of  these 
vibrations  In  the  molecule  will  shift  the  location  of  the  minimum 
density  of  states  towards  smaller  values  of  the  reaction  coordinate  and 
alternatively,  stiffening  of  Internal  modes  will  shift  the  minimum  density 
of  states  to  larger  values  of  the  reaction  coordinate  as  compared  to  the 
locution  of  the  centrifugal  barrier.^*'  As  the  temperature  goes  up,  the 
minimum  density  of  states,  as  well  as  the  top  of  the  rotational  barrier 
move  to  smaller  values  of  the  Inter  fragment  distance.  The  lumd  scission 

In  I ,5-hexadleno  to  two  resonance  stablll/.ed  allyl  radicals  Is  a special  case 
In  that  the  anticipated  loosening  of  Internal  modes  upon  bond  breaking  Is 
accompanied  by  n concomitant  "resonance  stiffening,"  mainly  of  the  free 
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lateriuil  rotutluns,  resulting  In  a near  cancellation  of  the  temperature 
ilepunilence  of  the  location  of  the  minimum  density  of  states.  The  near 
balance  of  two  opposlitg  effects,  namely  loosening  and  stiffening  of  Internal 
modes  In  the  critical  configuration  of  1 , 5-hexadlene  upon  bond  breaking, 
and  Its  apparent  lack  of  temperature  sensitivity,  resembles  the  temperature 
Independent  vibrational  model.  This  balance  Is  thus  thought  to  be  the  main 
reason  for  the  very  small  "tightening"  of  the  transition  state  with  Increasing 
temperature  as  reflected  in  the  Increase  of  H by  only  3%  from  31)0  to  900  K. 

The  e.xplanatlon  put  forward  above  Is  of  a qualitative  nature  and  warrants 
further  detailed  calculations  on  the  I ,5-hexadleno  system.  The  foregoing 
discussion  makes  It  clear  why  a vibrational  model  Is  ecjually  successful  In 
describing  the  bond  scission  of  I ,5-hexudlene . In  tltls  case,  both  the 
vibrational,  and  rotational  models  are  characterized  by  the  same  (negative) 
temperature  dependence  of  their  respective  Arrhenius  parameters,  though  to 
a different  dtjgree  (see  Table  V),  whereas  In  the  case  of  a C-C  bond  split  In 
an  alkane,  the  vibrational  model  predicts  a strong  positive  temperature 
dependence  ot  its  Arrhenius  parameters. 

A glance  at  the  high-pressure  Arrhenius  parameters  and  log  A^j  tor 
model  (H) , Table  V,  reveals  that  the  rotational  model  has  only  4 * tt/2  heat 
capacity  associated  with  Its  hindered  rotations,  such  that  and  log  A^ 
ilecreast-  by  2.7  leal  mol  amt  1.24  logarithmic  units,  respectively,  t>ver  tlie 
temperature  range  3l)l)  to  9t)l)  K.  Given  the  under lyliig  assumption,"^  tins 
IS  •.i»Mii«liat  in  illsaccoial  .^ith  ihe  result  from  the  t!qi»i  librium  sludt  (vide  supi.i), 
I’lch  predicts  a constan.  value  lor  - - dll")  over  this  lempeiMture  range 

Ith  the  present  method  ol  calculation  (l.e.,  choosing  an  A-factor  it  T = .100  K 
lO  match  the  low-temperature  recombination  rate  constant),  one  la  led  to  con- 
lode  that  the  system  has  a n€*gatlve  activation  energy  (or  recombination 
iK,.  * - 2.70  kcal  mol).  Instead  of  t 1.80  kcal  'mol  (RT)  as  required  by  the 
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This  Is  an  artifact  of  the  iiutdel  , howev<>r,  cind  is 


r 


HI  derlying  assumption.^’' 
clearly  the  consequence  of  lacking  heat  capacity  of  the  four  hindered  rotations. 

Ii  conversely,  one  proceeds  to  natch  the  high-temperature  fall-off  data  at  900  K 
with  = + 56.40  kcal/mol,  one  is  led  to  predict  negative  values  for  '1  at 
T = 300  K,  (E^  = 59.34  kcal/mol),  which  is  physically  unreasonable.  No 
attempt  has  been  made,  however,  to  associate  an  increased  amount  of  heat 
eapacliy  with  these  hindered  internal  rotations,  knowing  that  the  upper  limi* 
with  respect  to  the  heat  capacity  of  the  activated  complex  Is  represented  in 
the  vibrational  model  (1) . With  the  present  fall-off  data  for  recombination 
of  allyl  radical,  no  unambiguous  choice  between  transition  state  models  (1) 
or  (2)  can  be  made.  We  favor  at  this  point  the  simple  vibrational  transition 
state  model  (1),  because  it  doesn't  suffer  from  the  shortcomings  cited  above. 

Tsang**^  estimated  the  Arrhenius  parameters  for  the  bond  scission  of  BA  at 
1100  K (log  Aj  = 14.20,  E^  = 59.3  kcal/mole) , through  the  application  of  the 
usual  geometric  mean  rule  for  the  cross  combination-to-combination  ratio 
r(r  = kj.(AB) /(kj.(AA)kj.(BB) ) ).  Using  the  computed  value  of  14.20  for  log  Ajj , 
together  with  the  overall  entropy  change  (AS°)  of  —33.50  e.u.  at  T = 1000  K, 
yields  log  k^  = 8.79  (Table  VI,  reaction  4),  which  is  one-to-one-and-a-hal f 
orders  of  magnitude  too  low  in  light  of  our  results.  The  measured  recombination 
rate  constant  kj.  at  1000  K with  AS°  = —33.50  e.u.  for  equilibrium  (2)  results 
in  log  Ajj  of  15.75  or  15.26,  according  to  the  vibrational  or  rotational  transition 
state  model.”®  If  the  rate  constant  for  allyl  + CHg ‘ recombination  is  set  equal 
to  10*®’”°  (see  Table  VI,  reaction  2),  one  would  predict  a value  of  .95  or  1.68 
for  r,  instead  of  2.45  (as  used  by  Tsang)  on  the  basis  of  our  measured  value  for 

OD 

kj. , which  is  subject  to  the  choice  of  the  appropriate  transition  state  model. 

We  favor  a value  of  r = 0.95  as  a result  of  our  preferred  choice  of  the 
vibrational  model  (vide  supra).  A recent  pyrolysis  study  of  1,1  -azoi sobvitane” ” 
provided  an  experimental  value  of  r = .25  for  the  methyl  + isobutyl  radical  pairs. 
The  results  of  the  flash  photolysis  study®  of  2-methylbutene-l , where  r = .09 


EXPERIHENTAL  AND  CALCULATED  RECOMBINATIGM  RATE  CONSTANTS 
RADICAL  RECOMBINATION  REACTIONS  INVOLVING  STABILIZED  RADICALS 
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was  found,  seem  too  low  and  may  indicate  an  experimental  problem.  Adjust- 
ment of  r to  higher  values  in  this  case  would  decrease  kj.  for  reaction  (3) , 
Table  VI,  which  seems  somewhat  high  in  view  of  the  results  for  allyl  radical 
recombination  of  van  den  Bergh  and  Callear®  at  300  K,  and  concomitantly  would 
increase  for  reaction  (2),  Table  VI,  which  seems  unusually  low  for  the 
cross-combination  rate  constant.  These  few  examples  show  that  r = 2 .0  may  not 
always  be  valid  without  prior  experimental  assessment  of  the  involved  rate 
constants . 

We  conclude,  therefore,  that  the  present  "high"  value  for  the  recombin- 
ation rate  constant,  which  is  quite  similar  in  magnitude  to  common  alkyl 
radical  recombination  rate  constants  (ethyl  radical  recombination,  kj.  = 

(7.80  ±.  1.80)  * 10*  M~‘  s~* , ref.  33)  is  indicative  of  the  fact  that  the 
recombination  reaction  shows  no  apparent  activation  energy  and  that  the 
delocalization  of  the  unpaired  electron  has  no  effect  on  the  recombination 
rate  of  allyl  radical. 


V SUMMARY 


This  study  yields  (allyl)  = 38.3  ± 1.5  kcnl  mol~‘ , tind  thus 

I ,300 

an  allyl  resonance  energy  ARE  = 12.5  ±.  2.0  Real  mol~' . At  the  same 
time,  these  results  yield  a value  for  k“  which  Is  about  the  same  as  for 
alkyl  radicals,  indicating  no  effect  of  the  electron  delocalization  on 


reactivity . 
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CAPTIONS 


Figure  1 Block  illugriim  for  experimental  setup  of  V1J*P  (Very  Low-Pressure 
Pyrolysis)  experiment  . 

Figure  2 Very  Low-Pressure  Pyrolysis  Reactor  (133.5  cm^ ) 

Figure  3 Van't  Hoff  Plot  of  Equilibrium  Conslani  K = k,./k,j 

r , il  ‘ " 

k 


Figure  4 Brauman  Plot  (see  text)  for  the  Recombination  Reaction  k,. 

k_  ^ ^ 

r \ 

2 • ■ * a t 625  K 

(experiments  numbers  2 and  3,  Table  III) 


Figure  5 Experimental  and  Calculated  Rate  Constants  k,.  for  Recombination 
of  .Allyl  Radicals  In  Function  of  Temperature 

m Reference  8 
® This  work 
Vibrational  Model 

Rotational  Model  with  Hindrance  ( '!)  to  2-dimenslonal 

rotation  ns  Indicated  In  Table  V.  The  shaded  area 
corresponds  to  values  of  H ±.  5‘r 
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TAe  steady-state  kinetic  expressions  for  the  reaction  system  (1)  to  (5) 
In  a stirred  flow  reactor  shall  now  be  derived. 


Under  steady-state  conditions,  the  following  expressions 
(Al)  to  (A4)  for  the  concentration  of  the  Involved  species 
result ; 


(Al) 


d(DAO)/dt 


DAO 


ki(UAO) 


e 


(A2) 


d(CO^)/dt 


CO 

2k,  (DAO)  - k'"  *(CO,)h  0 
* e ^ 


(/.)) 


(A4) 


dCCjHjO/dt  = 2k,  (DAO)  + 2k(,(BA)  - 2k,.(C3il5') 

.C3H; 

‘e 

. BA 

d(ta)/dt  = kr(C3H,-)‘  - k^'AA)  “ k^  (BA)£0 


- k„*  ^(03H3-)5  0 


where  Rp  Is  the  flow  of  the  species  P Into  the  reactor  (In  units  of  molecules 

g-i  is  the  escape  rate  constant  of  P out  of  the  reactor,  and  DAO  stands  for 

e 

dlallyloxalate  and  BA  for  1 ,5-hexadlene . 

Rearrangement  of  (A4)  yields  (A5) ; 


(A5) 


BA 

(BA)/(C,H  •)“  = k^/(k^  + k^  ) 


Keeping  in  mind  that  R”  = . where  Rp  Is  the  flow  of  species  P out  of  the 

reactor,  the  combination  of  (A3)  and  (A4)  results  In  (A6)  using  the  substitution 


2k, (DAO)  = R°  from  (A2) ; 

* CO3 


(A6) 


(BA) 


(BA) 


’ CO,  BA  \2 

k * (CO3)  - 2k  (BA)  ) 
© © / 
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Comparison  of  (A6)  and  (A5)  results  In  the  desired  expression  (A7)  , which 
enables  one 


(A7) 


k 

r 


k + 
d 


3A 


2R° 

BA 


<=02 


2R“  V 
BA 


to  establish  k /(k,  + k„  ) with  the  experimental  values  of  R°  and  R“  . 

r (t  e BA  CO.. 

BA  ^ 

With  two  Independent  values  of  k (e.g.,  corresponding  to  a large  and  a 

6 

small  aperture  In  the  VLPP  reactor) , the  algebraic  separation  of  k^  and 
k^  Is  possible. 
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Chapter  2 


ABSOLUTE  RATE  CONSTANTS  FOR  METATHESIS  REACTIONS  OF 
ALLYL  AND  BENZYL  RADICALS  WITH  HI(l)I).  THE  HEAT  OF 
FORMATION  OF  ALLYL  AND  BENZYL  RADICALS^ 
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ABSTRACT 

^3 

The  metathesis  reaction,  Cglig * + HI(DI) ► CjHgCCgHgD)  + I,  has  been 

studied  In  the  gas  phase  using  the  VLPP  technique.  The  result  Is 
logCkg/M"^  s"^)  = (9.73  + .21)  - (4.0  ±.  1.0)/e  at  T = 1000  K and 

log(k3/M“^  s-^)  = (9.58  + .35)  - (3.0  ± l.O)/0  at  T = 635  K,  where  9 = 2.303 

RT  In  kcal  mol“^ . The  result  for  the  metathesis  reaction  CgHgCHg’  + DI  i 
CgHgCHjD  + I'  is  log(k3/M"*  s"*)  = (9.93  ± .22)  - (4.0  + 1.0) /0  at  T = 965  K. 

These  rate  expressions  were  extrapolated  to  lower  temperatures  using  a 

transition  state  model  In  order  to  compute  the  equilibrium  constants  for  the 
above  metathesis  reactions  using  the  rate  constants  for  the  reverse  metathesis 
from  classical  lodinatlon  studies.  The  equilibrium  constants  yield 
AH^(allyl)  = 39 .1  £ 1 .0  kcal/mol  and  AH^(benzyl)  = 46.60  ±1.5  kcal/mol  at 
T = 300  K.  These  values  correspond  to  stabilization  energies  of  11.7  ±1.5 
kcal/mol  and  11.3  + 2.0  kcal/mol,  respectively  (l.e.,  DH(allyl-H)  = 86.3  ± 1.0 


kcal/mol  and  DH(CeH5CH2-H)  =86.7+1.5  kcal/mol). 


I INTRODUCTION 


A large  number  of  values  for  free  radical  heats  of  formation 
(AH^  (R*))  have  been  obtained  by  the  spectrophotometric  iodination  technique 
which  is  well  documented  in  the  literature.^ 


The  pertinent  reactions  are: 


I2  + M 

21-  + M 

RI  + I' 

; ^12" 

R-  + HI 

ka 

RH  + !• 

; K3  ^ = *"3/k 

k. 

with  the  following  overall  equilibrium 


(1) 

(2) 

(3) 


RI  + HI  -r—  RH  + I2  (4) 

In  general,  the  temperature  dependent  values  of  kj  and  kg /kg  (starting 

with  RI  + HI)  or  and  ^ (starting  with  RH  + I^)  can  be  obtained.  The 

usual  method  of  extracting  AHj(R‘)  from  these  studies  is  to  assume  that 

E2  = 0.0+  1.0  kcal/mol  and/or  = 1.0  ± 1.0  kcal/mol.  (The  measured  values 

of  E2-E2  are  not  inconsistent  with  these  assumptions.)  Thus,  a measurement  of 

will  yield  a value  for  AH^  ^ and  AH^(R')  will  follow  since  the  appropriate 

2 

values  for  HI,  I*  and  RH  are  known. 


Very  low-pressure  pyrolysis  (VLPP)  allows  the  measurement  of  fast 
bimolecular  reactions  — such  as  (5)  in  the  gas  phase. 


R'  + HI 


RH  + I- 


(5) 


We  take  advantage  of  the  fact  that  the  measured  rate  constants  kg  are 
unencumbered  by  competitive  radical-radical  recombinations 

■ ji*. 


PiOB 
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relative  to  which  the  overwhelming  majority  of  fast  blmolecular  radlcal- 

3 

molecule  reactions  have  been  measured.  The  availability  of  Arrhenius  para- 
meters for  forward  and  back  reaction  rate  constants  provides  values  for 
AHj(R*)  and  entropies  of  radicals  (S°(R*)).  Accordingly,  we  report  in  this 
work  the  rates  of  reactions  (6)  and  (7) 


+ HI(DI) 


C-H.Cri-  + HI(DI) 


C^H^CC^HjD)  + V 


CsHsCHgCCeHsCH^D)  + l' 


together  with  their  Arrhenius  parameters  over  an  extended  temperature  range, 

thereby  obtaining  values  for  M°(R‘)  of  ally!  (CjHg’)  and  benzyl 

(C,H,CH-)  radicals  by  using  the  results  of  previously  published  determinations 
6 3 2 

4 5 

of  k, , A,  and  E,  from  lodlnation  studies  of  propylene  and  toluene. 

4 4 4 
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II  EXPERIMENTAL 


The  description  of  the  VLPP  molecular  beam  sampling  apparatus,  together 
with  the  all-quartz  reaction  vessel  (Knudsen  cell)  has  been  presented  in  a 
previous  publication.®  The  source  and  purification  of  the  radical  precursors 
dlallyloxalate  (CgH,j,0,),  3 ,3 ^-azo-l-propene  (CgHj^N^)  and  blbenzyl 
have  also  been  described.®  Benzylvlnylether  (CgHj^O)  was  synthesized 
according  to  literature  methods.^  GC-MS-analysls  showed  a 2%  impurity  of 
benzaldehyde  (C^HgO) , which  did  not  interfere  with  our  measurements.  The 
standard  gas-handling  system  was  heated  in  order  to  increase  the  vapor 
pressure  of  dlallyloxalate,  bibenzyl  and  ben  .ylvinyl  ether,  such  that  suitable 
flow  rates  of  the  radical  precursors  into  the  heated  Knudsen  cell  could  be  obtain- 
ed. Hydrogen  iodide  (Linde  Air  Products,  Inc.)  and  deuterium  iodide  (99  atom  %, 
Merck,  Sharp  and  Dohme,  Ltd.,  Canada)  were  used  without  further  purification. 

The  experiments  were  carried  out  by  setting  the  flow  of  the  radical  precursor 
to  a low  and  constant  value  and  monitoring  m/e  = 43  (CjH^D)  and  m/e  = 92  (C^Hg) 
or  m/e  = 93  (C^H^D) , as  a function  of  flow  rate  of  HI  or  DI . The  reaction  vessel 
had  two  inlet  (capillary  or  needle  valve)  systems  from  two  independent  gas-handling 
systems  so  that  the  components  met  for  the  first  time  in  the  hot  reactor. 


i 
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Ill  RESULTS  AND  DISCUSSION 


A.  Allyl  Radical  + HI (PI) 

Dlallyloxalate  (CgHjo04)  was  used  as  a precursor  for  allyl  radicals 
at  <T>  = 1000  K and  3 ,3  ^-azo-l-propene  (CgHj^N^)  at  <T>  = 750  K.  By 
using  two  different  precursors  for  allyl  radical  It  was  possible  to  cover  a 
range  of  380  K In  the  study  of  the  metathesis  (6) . Application  of  the  steady 
state  assumption  to  reaction  (6)  In  a low-pressure  stirred  flow  reactor  results 
In  the  following  relation®  (see  appendix  for  details  of  the  derivation) : 


1/f  = 1 + 


e e 

k,(HI) 


(8) 


where  f Is  the  fraction  of  radicals  ’’titrated"  at  a certain  steady-state 

, X C,H. • HI  ^ , 

concentration  (HI)  , k , and  k are  the  escape  rate  constants  (s”  ) 

e e 

of  allyl  radical  and  HI  out  of  the  reactor,  f Is  defined  as  (C3Hg)/(C3Hg)^, 
where  (C3Hg)  Is  the  steady-state  concentration  of  C3Hg  at  a certain  (HI),  and 
(C3Hg)^  Is  the  same  concentration  at  (HI)  = where  essentially  all  radicals 
have  reacted  with  HI . In  the  case  of  dlallyloxalate  as  precursor  for  allyl 
radical, (CO3)  could  be  taken  as  representative  of  the  amount  of  the  allyl 
radicals  present  In  the  reaction  system,  so  that  1/f  could  be  equated  to 
(CO  )/(C3Hg)  or  R°  . This  method  proved  to  be  valid  because 

COj  CjHg 

plots  of  1/f  vs.  R°  yielded  straight  lines  with  unit  Intercepts,  confirm- 

COj 

Ing  that  the  amount  of  CO,  present  was  Indeed  representative  of  the  amount  of 
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allyl  radical  In  the  system.  The  following  assumptions  have  been  made  In 

deriving  equation  (8);  (HI)  » (CjHg*)  and  = F°  , where  F°  is  the 

flow  of  HI  in,  and  F°^  is  the  flow  out  of  the  reactor.  Under  our  experimental 

conditions  both  assumptions  were  found  to  be  Justified,  so  that  a plot  of  1/f 

vs.  1/(HI)  gave  a straigSttt  line  with  unit  intercept  (within  experimental  error) 

C H * HI 

and  )/h3  as  slope,  from  which  kj  was  readily  obtained  given  the 

reactor  parameters.^  Special  attention  was  paid  to  the  requirement  of  carrying 
out  the  titration  at  the  lowest  possible  concentration  of  C3H3 * in  order  to 
suppress  the  bimolecular  recombination  reaction.  Figure  1 displays  a typical 
plot  1/f  vs.  1/(HI)  and  Table  I lists  the  rate  constants  k3/M“‘  s“‘  as  a 
function  of  temperature . These  are  plotted  in  the  usual  Arrhenius  form  in 
Figure  2.  The  results  are: 

kg  = 7.2  X 10®  M~‘  s”'  at  <T>  = 1000  K 

k3  = 3.6  X 10®  M~‘  s~‘  at  T = 635  K 

The  Arrhenius  plot  of  Figure  2 accommodates  slight  curvature  over  the  temperature 

range  635-1000  K.  For  very  exothermic  reactions,  such  as  (6)  and  (7),  with 

very  small  activation  energies,  one  expects  a small  isotope  effect  with 
HD  _ 

kj  /k3  = 1.5  ±.  .5.®  This  ratio  of  the  rate  constants  can  be  seen  at  higher 

temperatures  (Figure  2),  where  both  HI  and  DI  were  used  as  titrating  agents. 

It  may  be  seen  that  E3  at  <T>  = 1000  K is  4.0  + 1.0  kcal/mol  (log  A3/^“*  s“^  = 
9.58  + .35),  significantly  different  from  the  "usual"  assumed  value  of  1 .0  i 1 .0 
kcal/mol  used  throughout  the  literature.* 


Table  1 


RATE  C(»STANTS  IXIR  METATHESIS  REACTIONS 
AS  A KUNCTK*  OK  TEMPERATURE 


tLk] 

kjLM-*  s-‘ 

k., 

V + DI  i — 1 

. I- 

907 

4.80  * 10 

911 

4.90  * 10 

929 

6.68  * 10' 

929 

6.99  * 10' 

926 

6.37  10' 

1007 

8.06  » 10' 

10 1-1 

6.76  10' 

. HI  — i*. 

635 

3.57  * 10' 

697 

5.50  * 10' 

761 

4.83  * 10' 

878 

7.24  * 10' 

982 

8.54  10' 

a 

a 

a 

a 

a 

a 

a 


b 

b 

b 

b 

b 


**Radical  source:  Diallyloxalate  CgHj^O^ 

b / 

Radical  source:  3,3  -azo-l-propene 
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et  al.  £o^-d  the  following  rate  expression  at  <T>  - 530  K: 
log(k^AJ-‘  s'M  = 10.25  - 18.04/d,  where  9 = 2.303  RT  in  kcal  mol"*.  in  order 
to  provide  a value  at  1000  K for  the  equilibrium  constant  K.,  .3  = *<4 /•‘j  13), 
where  R*  = CjHj",  a transition  state  model,  which  has  been  shown  by 
Benson  and  coworkers  to  be  appropriate^  in  cases  of  metathesis  reactions,  such 
as  (6)  and  (7)  , was  used  to  e.xtrapolate  the  measured  rate  constant  k.,  from  530  K 
to  1000  K.  The  transition  state  (9) 


for  reaction  (6)  or  Its  reverse  was  approximated  by  taking  allyl  iodide 

(C3H5I)  as  a model,  adding  two  CHI  bending  and  one  (CH)-I  stretching  mode  to 

CjHjI  and  including  other  minor  corrections  due  to  the  presence  of  the  extra 

H-atom  in  (9).  Table  II  demonstrates  the  method  and  summarizes  the  result 

for  the  entropy  of  activation  AS^’  ^’f  the  reverse  of  reaction  (6).  With 

= +3.94  e.u.,  the  following  rate  expression  at  T = 1000  K results: 
P.+ 

log(k^/M"'  s-‘)  = 11.32  - 19.90/6,  T = 1000  K. 

K.  , at  T = 1000  K then  turns  out  to  be  1.39  * lO"'*  which,  together  with 
an  estimate  for  S^lallyl  ) and  Cp  data''*’*”  yields  a "third  law"  standard 
heat  of  formation  AH^lallyl  ) *39.1  +1.0  kcal/mol  at  T = 300  K.  The 

'second  lav"  heat  of  formation  derived  from  the  comblna t ion  of  the  .Arrhenius 

^The  combined  error  limits  have  been  estimated  to  be  ±.  1.0  kcal  mol  in  view  tf 
the  high  precision  determinations  of  k^  Iref.  4)  and  k3  (this  work). 


52 


Table  rr 


ESTIMATED  ENTROPY  OK  ACTIVATION  FOR  TKE 
METATHESIS  REACTION,  CjHf,  + l'  C3H5  • + 

(Standard  State:  1 atm) 


Corrections 


30  ^ - c®  r ■ 


s°cc-H,,i)  - s^cr ) 


- S"(C3He) 


300  K 

530  K 

1000  K 

-29.31 

-30  .05 

-31  .17 

Corrections : 

spin  (Rln2) 

1 .39 

1 .39 

1 .39 

rotation^ 

H 

1 .69 

1 .69 

1 .69 

2 bending  (C"  I):  300  cm~‘ 

2 * 1 .40 

2 2.40 

2 * 3.60 

H 

1 stretch  (C  I) tlOOO  cm~* 

.10 

.50 

1 .40 

resonance  stiffening  of  internal 

- 1.30 

- 1.43 

- 1.12 

rotation  = 4 -•  13  kcal/mol 

C=C,  1650  cm-’  C-C,  1400  cm"’ 

0.0 

0.0 

+ .10 

C-C,  1000  cm-’  - C-C,  1400  cm-’ 

- .10 

- .30 

- .50 

420  cm-’  -♦  635  cm-’ 

- .52 

- 1.00 

- .90 

-25  .52 

-24.40 

-21  .90 

Adjusted  to  give  experimental  at  <T>  = 530  K (log  A^/M”*  = 10.25) 


expressions  is  In  excellent  agreement,  since  the  value  of  S^iallyl)  which 
is  ilerived  from  A3  and  is  62.80  e.u.,  compared  to  the  value*”  of  62.10  e.u. 
used  in  the  "third  law"  calculations.  Table  111  summarizes  the  kinetic  and 
thermochemical  results  and  provides  a good  example  of  an  Internally  consistent 
set  of  Arrhenius  and  thermochemical  parameters. 

The  allyl  radical  heat  of  formation  derived  above  from  two  separate  exper- 
iments is  in  excellent  agreement  with  a recent  equilibrium  study  and  other 
determinations  of  AHjCallyl  ) . ” The  present  value  for  iiaijCallyl  ) provides 
an  allyl  resonance  stabilization  energy  (ARE)  of  11.7  i 1 .5  kcal/mol  (see,  e.g., 
reference  11  for  the  definition  of  ARE).  The  large  variation  of  Aj  and  A^  with 
temperature  (Tabic  III)  is  noteworthy;  this  brings  about  a concomitant  change 
in  E3  and  E^  which  tend  to  partially  cancel  the  increase  In  the  A-factor. 
Obviously,  this  cancellation  of  log  A versus  E is  only  partial,  so  that  the 
suggested  curvature  in  the  Arrhenius  plot  in  Figure  2 results.  It  has  repeatedly 
been  pointed  out  in  the  literature,*  that  the  Arrhenius  A-factor  tor  the  reaction 
of  1’  with  an  alkane  or  olefin  Is  quite  high  (log  A4  = 10.3  at  550  K) . Such  a 
behavior  is  obvious  from  Table  III  in  the  present  case  as  well. 

n . Benzyl  Radical  + HKUl) 

At  <T>  = 880  K benzyl  radicals  were  generated  by  unlmolecular  decom- 
position of  ben  y Ivlny lether  (C^H^-OC^Hj)  and  at  <T>  = 1060  K by  the  bond 
breaking  process  of  bibenzyl  (C^H,-C,H,).  The  evaluation  of  the  rate  constant 
k3  (7)  follows  the  same  linos  ns  in  the  case  of  allyl  radical  . Figure  3 displays 
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Table  III 


a typical  1/f  vs.  1/(HI)  plot  (equation  8),  where  f represents  by  analogy 
with  the  definition  given  In  section  A,  the  fraction  of  radicals  titrated  at 
a certain  steady-state  concentration  of  HI  with  respect  to  the  total  concen- 
tration of  radicals  present  in  the  absence  of  HI;  Table  IV  lists  the 
resulting  rate  constants  kj/M"*  s"*  as  a function  of  temperature,  and  Figure 
4 shows  the  corresponding  Arrhenius  plot.  By  using  two  different  benz;yl 
radical  precursors,  which  decompose  at  different  temperatures,  it  was  possible 
to  study  the  metathesis  reaction  (7)  over  a range  of  315  K. 

The  rate  constant  kj/M"*  s“'  was  found  to  be  1.08  n 10®  M”'  s”*  at 
<T>  = 965  K.  The  activation  energy  Eg  at  965  K is  4.0  ±.  1.0  kcal/mol 
(log  Ag/M"*  s~‘  = 9.93  + .22)  in  remarkable  agreement  with  the  parameters  found 
for  allyl  radical.  The  Arrhenius  plot  in  Figure  4 suggests,  by  analogy  with  the 
allyl  case,  slight  curvature  at  higher  temperaiures  and  a possible  small  sys- 
lematic  error  between  the  two  series  of  different  precursors  for  benzyl  radical, 
but  it  is  obvious  as  well  that  the  straight  line,  corresponding  to  Eg  = 4 .0  ± 1 .0 
kcal/mol  does  not  do  violence  to  the  present  data.  Due  to  the  more  preliminary 
nature  of  the  iodlnation  study  of  toluene®  the  experimental  rate  of  the  reverse 
of  reaction  (7)  at  T = 500  K was  assumed  to  be  correct  (1.41  10“  M"‘  s“‘ ) 

rather  than  using  the  published  Arrhenius  parameters,  since  A^/M"*  s"‘  seems 
to  be  too  low  by  at  least  an  order  of  magnitude.  In  order  to  derive  a value 
for  the  equilibrium  constant  g = k^/kg  at  T = 500  K (3),  where  R*  = C.,H^  ’ , 
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Table  IV 

RATE  CONSTANTS  1X)R  THE  METATHESIS  REACTION 
FUNCTION  OF  TEMPERATURE 


k,TM-‘  S-* 


7,08  10® 

5.78  * 10® 
9.03  « 10® 
6.99  10“ 

9.35  » 10® 

1.35  * 10® 


1,23  * 10® 
1 .40  * 10® 


1.51  10® 
1,27  * 10® 
1 .70  * 10® 
1.37  * 10® 
1.32  10® 
1.72  10® 
2.30  * 10® 


Radical  Source:  Benzylvlnylether  (C^HgCH^OC^Hj ) 
Radical  Source:  Blbenzyl (1 ,2-Dlphenylethane) 


1 ••  wr 


the  transition  state  model  (10)  was  used  by  analogy  with  the  allyl  case: 


(10) 


discussed  above.  Benzyliodide  (C.,H,1)  was  taken  as  a model  and  corrected 
for  the  transition  state  (10).  Table  V shows  the  detailed  corrections 
and  the  computed  entropy  of  activation  for  reaction  (7).  The  rate 

expression  for  kj  extrapolated  to  T = 500  K is: 

logikj/M-*  s-‘)  = 9.00  - 2.H1/0  Ul) 

3 at  T = 500  K is  then  computed  to  be  2.34  x 10"*’,  which  yields  a 

standard  heat  of  formation  ^j(benzyl)  = 46.60  t 1.50  kcal/mol  at  300  K, 

using  an  estimate  for  S°(benzyl)*“  and  Cp  data.*'*” 

The  comparative  rate  shock-tube  decomposition  of  i sobutylbenzene* * at 
<T>  = 1100  K yields  AHj(benzyl)  = 48.42  ±2.0  kcal/mol  at  300  K under  the 
assumption  of  zero  activation  energy  (in  pressure  units)  (or  the  reverse 
reaction  at  reaction  temperature  ( Arrhenius  activation  enei’gy  for  decom- 
position = ah”),  a value  markedly  higher  than  the  currently  accepted  value 

1 ^ 

from  lodination  studli?s  and  related  evidence  of  45.00  ±1.5  kcal  mol.  "“lie 


b 

Assumption  of  zero  activation  energy  in  concentration  units  at  <T  > = 1100  K 
would  result  in  an  even  higher  value  (or  AH”(benzyl)  at  T = 300  K (50.60 
kcal/mol) . 
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Table  V 

ESTIMATED  ENTROPY  OF  ACTIVATION  TOR  THE 
METATHESIS  REACTION,  C^H^  • + DI  C^H^D  + I 
(Standard  State:  I atm) 


CH^-  r CH^  D I]  CHgl 

ip'  + DI  0 - (Q)  + 

As"'*'  = S^CC^H^I)  -S"(DI)  - S^CC^H,  ■) 


= LUJ  + corrections 


±' 

Asp 

Corrections : 

spin  (Rln2) 

rotation® 

D 

2 bending  (C 

I)  ; 300  cm”‘ 

1 stretch  (C 

I)  : 1000  cm"‘ 

Resonance  stiffening  of  Internal 
rotation,  = 3 -•  15.5  kcal/mol 

C -C,  1000  cm-‘ 

- C - C , 1200  cm-* 

ASp 

300  K ! 500  K 965  K 


- 34.95  I - 36.32  | - 37.89 

I I I 

i ' i 

1 .39  I 1 .39  j 1 .39 

1 .98  ! 1 .98  1 .98 

) ' 

2 * 1.40  i 2 * 2.30  2 * 3.53 

I 

.10  J .5 

- 2.10  I - 1 .0 


II 


- 30.86 


Adjusted  to  give  experimental  A3  at  <T>  = 965  K (log  A3/M 


discrepancy  of  the  results  of  both  experimental  methods  is  not  clear.  Our 
value  for  AHj(benzyl)  lies  right  in  between  the  two  seemingly  different 


results,  and  both  seem  compatible  with  our  value  from  this  study  given  our 
overall  experimental  error  of  + 1 .5  kcal/mol . This  result  is  quite 


satisfactory,  given  that  only  the  rate  for  the  reverse  of  (7)  was  used.  A 
consistency  check  on  the  standard  entropy  of  benzyl  radical  using  the 
Arrhenius  parameters  Aj  and  A^  (2.30  R log  A4/A3  = AS® ^3)  could  not  be 
performed  because  of  the  lack  of  reliable  Arrhenius  parameters  for  the 
reverse  of  (7)  Instead,  Arrhenius  parameters  were  calculated  for  the 
reverse  of  (7)  using  an  estimate  for  AS^^a*’'®  and  the  experimental  A3 
together  with  the  rate  at  T = 500  K.  These  results  are  listed  in  Table  VI. 
The  same  observations  can  be  made  as  in  the  allyl  case  when  considering  the 


magnitude  of  log  A4  and  log  A3 , respectively. 


Table  VI 


MEASTOED  AND  CALCULATED  RATE  PARAMETERS  AND  OVERALL 
ENTROPY  CHANGE  FOjR  THE  EQUILIBRIUM* 


Toluene  + I ^ benzyl  + HI 


T/®K  ! logCkg/M'*  s-^)  I log(k4/M-^  s"^) 


AS^  3/e.u. 


’ E4  obtained  through  ACp  ^ = 3.60  e.u.  and  E4  at  T = 500  K (determined 
from  experimental  rate  of  reverse  of  (7)  and  A4) ; A4  calculated  from 
A4  at  500  K and  ACp  ^ . 

A4  from  A3  and  ASj  3 . E4  determined  from  experimental  rate  of  reverse 
of  (7)  and  A4 . 


Underlined  values  represent  experimental  results;  others  are  calculated 
Arrhenius  parameters  (using  transition-state  model  (10)  and  Table  V , 

0 = 2.303  RT  in  kcal/mol  ) . 


i 


60 


IV  KKFERENCES 


1.  D,  M.  Golden  and  S.  W.  Benson,  Chem.  Rev.,  125  (1969). 

2.  S.  W.  Benson,  Thermochemical  Kinetics,  2nd  Ed.,  John  Wiley  and  Sons,  Inc., 

New  York,  1976. 

3.  J.  A.  Kerr  in  Free  Radicals,  Vol . 1,  J.  K.  Kochi,  Editor,  John  Wiley  and 
Sons,  Inc.,  New  York,  1973 

4.  D.  M.  Golden,  A.  S.  Rodgers,  and  S.  W.  Benson,  J.  Amer.  Chem.  Soc . , 88 , 

3196  (1966)  . 

5.  R.  Walsh,  D.  M.  Golden,  and  S.  W.  Benson,  J.  Amer.  Chem.  Soc.,  650  (1966). 

6.  M.  Rossi,  K.  D.  King,  and  D.  M.  Golden,  preceeding  paper. 

7.  W.  H.  Watanabe  and  L.  E.  Coulson,  J.  Amer.  Chem.  Soc.,  79,  2828  (1956); 

H.  Yuke,  K.  Hatada,  K.  Nagata,  and  K.  Kajlyama,  Bull.  Chem.  Soc.  (Japan), 

3546  (1969). 

8.  N.  A.  Gac,  D.  M.  Golden,  and  S.  W.  Benson,  J.  Amer.  Chem.  Soc.,  9^,  3091  (1969). 


9.  A.  F.  Trotman -Dickenson , G.  S.  Milne,  Tables  of  Bimolecular  Reactions, 
NSRDS-NBS9,  National  Bureau  of  Standards,  Washington,  D.C.,  1967. 


10  . 

H. 

E.  O'Neal 

and 

S . W . Benson , Int . 

J.  Chem.  Kinetics,  1, 

217 

(1969) . 

11 . 

D. 

M.  Golden, 

N . 

A . Gac , and  S . W . 

Benson,  J.  Amer. 

Chem . 

Soc 

. , 91^,  2136  (1969) 

12  . 

W . 

Tsang,  Ini 

. J 

. Chem.  Kinetics,  2 

, 311  (1970)  ; VV  . 

T-iang 

. In 

! . J . Chem . 

Kinetics,  10,  41  (1978). 


61 


L 


CAPTIONS 

Figure  1 Experimental  determination  of  kg  (see  equation  8)  for  the 

k-  . 

metathesis  reaction  allyl'  + DI  propylene— dj  + I at 

1014  K.  Representative  example  from  Table  I. 


Figure  2 Arrhenius  plot  of  rate  constants  from  Table  I, 

metathesis  reaction  involving  Dl 

metathesis  reaction  Involving  HI 


Figure  3 Experimental  determination  of  kg  (see  equation  A-10)  for 

. k. 

the  metathesis  reaction  benzyl  + DI  w toluene-dj  + I 
at  989  K.  Representative  example  from  Table  IV. 


Figure  4 Arrhenius  plot  of  rate  constants  from  Table  IV. 

■ bibenzyl  as  precursor 

• benzylvinylether  as  precursor  for  benzyl  radicals 


fRKSDUiO  Fid  lUHC 


ob/ 


2.303  RT 

3000  ± 1Q( 
2.303  RT 


I 


2.00 


F°o,'’  /(molec  s*’)"’  [lO-’^l 

TA-32252:-'9 
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Appuntllx 


For  a radical  titration  (metathesis)  with  HX,  where  X = 1,  Br,  the  per- 
tinent reaction  system  Is: 


(A-1) 


(A -2) 


(A -3) 


ir  + HX  ► UH  + X' 

t.K’ 


Under  steady  state  conditions,  the  following  expressions  (A— 1)  to  (A-6) 
for  the  concentration  of  the  Involved  species  result,  where 


(A-A) 


d(RH)/dt  = k„(ll*)(HX)  - (RH) 


(A-5)  d(R-)/dt  = Ri^  - kg  (R*)  - k./R* ) (HX)-O 

I HX 

(A-6)  d(HX)/dt  = Rj^x  “ ” k„(R')(HX)  O 

p 4 

kg  Is  the  escape  rate  constant  of  species  P and  Rp  Is  the  flow  rate  of  P Into 
the  VLPP  reactor  (In  units  of  molecules  s“‘  -t“‘).  With  the  use  of  (A-4)  to 

(A-6)  , (RH)  can  be  expressed  ns  a function  of  (HX) : 

k k R^ . R^ 

a a V R HX 

(-7)  (RH)  = — (R-)(HX)  = * — — — 

RH  RR  .R  . /..vrN  . R^  . 

k k k + k (HX)  k + k (R  ) 

e e e a e a 

(A-7)  Is  simplified  to  give  (A-8)  under  the  condition:  (HX)  » (R*): 


(A-8) 


(RH)„  = 
00 


k (HX) 
a 


fRICXDlMO  PiOB  ht  AKy 


Crh) 

with  the  definition  1/f  = , the  following  expression  (A-9)  is 


(RH) 


HX. 


r\  ii/v. 

obtained,  where  (HX) : 

HX 


(A-9) 


1/f  = 1 + 


R"  HX 
k k 
e e 

k R° 
a HX 


R° 

HX 


In  the  limit  of  high  (HX) , where  R^  - R° 

HX  HA 


holds  to  a good  approximation 


and  where  R^  is  negligible  with  respect  to  R°  and  R^^ , respectively;  (A-9) 
R HX  HX 

can  be  simplified  to  expression  (A-10)  , which  relates  the  experimental 

quantities  1/f  and  R^jj  to  the  desired  rate  constant  k^  if  the  escape  rate 
R * HX 

constants  k and  k are  known. 


(A-10) 


1/f  = 1 + 


R’  HX 
k k 
e e 

^HX 


70 


LJL 


Chapter  3 
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ABSTRACT 


The  thermal  unimolecular  decomposition  of  three  vinylethers  has  been 

studied  in  a VLPP  apparatus.  The  high-pressure  rate  constant  for  the 

retro-ene  reaction  of  ethylvinylether  was  found  to  be  log  (k/s“‘)  = 

43 .4  ± 1 .0 

(11.47  + .25)  — at  <T>  = 900  K and  that  of  t-butylvinylether 

2.303  RT 

. 38  .4  + 1 .0 

proved  to  be  log  (k/s“*)  = (12.00  ± .27)  — — ■ at  <T>  = 800  K.  No 

^ .oDo  RT 

evidence  for  the  competition  of  the  higher  energy  homolytic  bond-fission 

process  could  be  obtained  from  the  experimental  data.  The  rate  constant 

obtained  for  the  C-0  bond  scission  reaction  in  the  case  of  benzylvinylether 
. 49  .3  + 1 .0 

was  log  (k/s-^)  = (15.05  t .30)  - o at  <T>  = 750  K.  Together 

2 .303  RT 

with  AH°  (benzyl")  = 46.6  kcal/mol,  the  activation  energy  for  this 

t »3  00 

reaction  results  in  M?  (CH.CHO)  = +3.9  ±.2.0  kcal/mol  and  in  a resonance 
stabilization  energy  of  3.2  + 2.0  kcal/mol  for  2-ethanalyl  radical. 
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I INTRODUCTION 


Most  vinyl  ethers  have  been  reported  to  decompose  by  molecular  elimin- 
ation mechanisms,  and  cyclic  transition  states  have  been  frequently  suggested 
as  a result  of  the  negative  entropies  of  activation  and  low  activation  energies 
observed.^  Recently,  it  was  observed  that  ethylvinylether  yielded  not  only 
acetaldehyde  and  ethylene  (c.f.,  (1)  ) , the  products  of  the  retro-ene 
molecular  elimination,^  but  also  significant  amounts  of  ketene,  ethane,  and 
butane  upon  multiphoton  dissociation  by  an  intense  CO^ -laser. ^ A higher 
activation  energy  reaction  path,  the  bond  sicission  process  (2)  , 


CgHgOCgHa  — l^CHgCHO  + 

(1) 

C2H5OC2H3  C^Hg • + CH^CHO 

(2) 

is  held  responsible  for  the  occurrence  of  the  latter  products,  which  are 
believed  to  be  due  to  secondary  reactions  of  radicals  produced  in  reaction 
(2) . Using  this  assumption,  an  apparent  reaction  "temperature"  of  ~ 1600  K 
could  be  inferred  from  the  branching  ratio  kj/kj , when  irradiated  with  the 
COg  laser. 

Specifically,  we  embarked  on  a complementary  thermal  study  of  the 
decomposition  of  various  vinyl  ethers  in  order  to  detect  and  study  the 
reactivity  of  the  radicals  C^Hg'  and  CHgCHO  (ethanal-2-yl)  with  respect  to 
ketene,  ethane,  and  butane  formation.  In  addition,  the  quantitative 
assessment  of  the  Arrhenius  parameters  for  the  bond  scission  of  benzylvinylether 
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to  benzyl  and  ethanal-2-yl  radical  at  high  temperatures  gives  the  heat  of 


formation  (AH$  ) and  resonance  stabilization  energy  of  ethanal-2-yl . We 

X|300 

report  the  unlmolecular  decomposition  of  ethylvlnylether  (CgHg-O-CgHg ) , 
t-butylvlnylether  (t-C^Hg-O-Cj^s)  i and  benzylvlnylether  (CgHgCHgOC^Hg ) In 
a VLPP  apparatus  up  to  temperatures  of  1300  K. 

II  EXPERIMENTAL 

The  molecular  beam  sampling  VLPP  apparatus  with  phase-sensitive  detec- 
tion has  been  described  In  detail  elsewhere,®  and  the  methods  of  assessing 
the  unlmolecular  reaction  rate  constant  from  the  mass  spectroscopic 

data  have  been  described  In  a review  article.^  The  VLPP  reactor  parameters 
are:  V = 133.4  ml , uj  = 4982  * (T/M)'/®  s"' , kgCB)  = 2.5571  * (T/M)‘/2  s-‘, 
1^(S)  = .2088  * (T/M)^/®  s"‘ . Ethylvlnylether  (Aldrich)  and  t-butylvlnyleth 
(Haven  Chemicals)  were  used  after  trap-to-trap  distillation.  GC-MS  analysis 
of  t-butylvlnylether  revealed  negligible  amounts  of  Isobutene,  t-butanol, 
and  t-butyl  acetate.  Benzylvlnylether  was  synthesized  following  standard 
procedures.®  (JC-MS  analysis  revealed  the  presence  of  ~ 2%  benzaldehyde , 
which  was  removed  by  distillation  in  vacuo.  The  low  vapor  pressure  of 
benzylvlnylether  made  necessary  the  use  of  a heated  inlet  system  when  higher 


flow  rates  were  required.  In  all  experiments  with  the  vinylethers,  except 


for  the  case  where  the  heated  inlet  system  had  to  be  used,  an  intei-nal  standard 


tCO^ , CHFj)  was  used.  Ueuterlum  Iodide  (99atom%)  was  purchased  from  Merck, 
Sharp,  and  Uohme  and  was  used  without  further  purification.  The  gas  mixtures 
were  prepared  in  a standard  gas-handling  vacuum  line  and  stored  in  S-f- 
darkened  bulbs. 

Ill  RESULTS 

A . Ethylvinylether 

For  the  decomposition  of  ethylvinylether,  the  pertinent  reactions  are 
(1)  and  (2).  The  course  of  reaction  (1)  was  followed  by  monitoring  the 
first-order  disappearance  of  ethylvinylether  (me  = 14,  15,  26,  27,  28,  29, 
31,  42,  43,  44,  45,  72)  at  m/e  = 72.  The  appearance  of  the  reaction 
products  acetaldehyde  (m/e  = 44,  43,  29)  and  ethylene  (m/e  = 28,27,26)  couM 
not  be  observed  ciuanti tatlvely , because  the  mass  spectrum  of  ethylvinylether 
had  a significant  contribution  to  the  intensities  of  almost  all  of  the 
peaks  of  the  reaction  products.  Ketene  was  not  detected  up  to  temperatures 
of  1270  K,  although  an  RRK  calculation  with  the  Arrhenius  parameters*' 
log  (k/s"*)  = 11.50  — 44.5/*),  predicted  detectable  yields  under  our  con- 
ditions (k,/kj,  ^ 5.5  at  T = 1300  K)  . Control  experiments  with  acetic 
anhydride  (3)  at  1270  K confirmed  that  ketene  (m/e  = 42,  41,14)  is 

CH3C0(X:0CH3  - CH3CO  + CH3COOH  (3) 

* , 

= 2.303  RT;  8 = 22,  u)  = 4982  x (T  Al)  ‘ ^ s"*  (reference  7). 
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stable  under  our  conditions,  and  also  produced  acetic  acid  (m/e  =60,  45,  15), 
as  well  as  methane  and  CO^ , which  apparently  stem  from  the  decomposition 
of  acetic  acid  on  the  walls  of  the  reactor.  It  is  noteworthy  that  acet- 
aldehyde decomposes  at  high  temperatures  to  give  CO  and  CH^ , probably 
according  to  (4) , since  the  presence  of  methyl 

CH3CHO  CHg'  + CHO'  - CH4  + CO  (4) 

radicals  was  established  through  titration  with  DI  to  give  easily  detected 
CHgD  (m/e  = 17,  16,  15).  The  absence  of  ketene  forces  us  to  conclude  that 
the  bond  fission  reaction  (2)  is  not  competing  under  our  experimental 
conditions  up  to  1270  K. 

The  rate  constant  kj  was  independent  of  the  flow  rate  over  the  range 

10^^  < S.  5 * 10*®  molecules  s"’ , and  the  experimental  values  of  k.,„</s“* 

EVE 

versus  T/K  are  plotted  in  Figure  1.  Moreover,  no  dependence  of  on  the 

aperture  size  of  the  VLPP  reactor  was  found,  further  confirming  the  absence  of 
secondary  reactions.  RRK  calculations  indicated  that  kj  is  in  the  fall-off  range 
in  the  temperature  range  studied  (750-1050  K) . In  order  to  obtain  the  high- 
pressure  Arrhenius  parameters  from  the  fall-off  data  of  Figure  1,  RRKM  calcu- 
lations were  performed  using  a fixed  transition  state  model  (see  Appendix  and 
Table  I)  and  varying  the  activation  energy  to  match  the  data.  The  previously 
measured  Arrhenius  A-factor*  log  (Aj/s~*)  =11.5  at  T = 800  K for  the 
six-membered  transition  state  of  the  retro-ene  reaction  (1)  was  assumed 
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Table  I 


MOLECULAR  PARAMETERS  ENTERING  THE  RRKM-CALCULATION  OF 
THE  RETRO-ENE  REACTION  OF  ETHYLVINYL ETHER.  C^HjOC^Ha  - C^^H^  + CH3CHO 


Molecule 

0 

Activated 

Complex 

1 1 

H 

3080(3) 

3080(3) 

2910(5) 

2910(5) 

1640(1) 

1640(1) 

1430(3) 

1430(3) 

Frequencies 

1320(3) 

1320(3) 

and 

1170(1) 

1170(1) 

1070(4) 

1070(3) 

degeneracies 

960(1) 

1050(1) 

800(3) 

960(1) 

580(1) 

800(2) 

350(4) 

700(1) 

200(3) 

600(2) 

100(1) 

400(2) 

350(4) 

10*“°I^I^I^  / (gr  cm*)* 

3,93  * 10® 

3.93  * 10® 

E°  / kcal/mol 

44.1 

®300  / ® -ii* 

80.25 

72  .83 

/ s-‘) 

11  .47 

Egoo  / kcal/mol 

43.4 

1o8<A3oo  / s-‘) 

11  .61 

CO 

E300  / Rcal/mol 

43  .7 

to  be  correct . 


Figure  1 displays  the  results  of  the  RRKM  calculation  (for 


reaction  1)  whose  high-pressure  Arrhenius  parameters  are: 

log  (ki/s-‘)  = 11.5  - 43.4/0  <T>  = 900  K (5) 

The  activation  energy  is  believed  to  match  the  data  to  within  + 1 kcal/mol 
with  a concomitant  uncertainty  in  log  (Aj/s“*)  of  ± .25. 

B . t-Butylvinylether 

In  order  to  observe  the  bond  sicission  reaction,  another  vinyl  ether 
with  a lower  activation  energy  than  EVE  was  chosen.  Thus,  t-butylvinylether 
was  decomposed  in  a VLPP  reactor.  The  pertinent  reactions  are  (6)  and  (7); 

t-C4H90CjjH3  -L^CHgCHO  + i-C^Hg  (6) 

t-C^HgOCjjHa  iil^CH^CHO  + t-C^Hg*  (7) 

The  extent  of  the  unimolecular  decomposition  was  observed  by  monitoring  the 
disappearance  of  t-butylvinylether  at  m/e  = 100  or  57  for  reasons  elucidated 
above.  Again,  no  indication  of  the  appearance  of  ketene  which  would  be 
indicative  of  the  occurrence  of  the  higher  energy  reaction  path  (7) , at 
m/e  = 14  or  42  could  be  Cound  up  to  temperatures  of  1220  K,  Indicating  that 
the  bond  scission  reaction  is  too  slow  to  compete  with  the  retro-ene  reaction 
Flow  rate  studies  of  the  unimolecular  reaction  rate  constant  k,  revealed 
no  dependence  in  the  range  8 .5  * 10*®  ^ ^EVE  ^ * 10‘®  molecules  s“* . The 

values  obtained  for  kyjjj  are  plotted  versus  temperature  in  Figure  2.  At  the 
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\ 

] 

temperatures  of  decomposition  (625-925  K)  , is  in  the  fall-off  range, 

so  that  the  respective  high-temperature  Arrhenius  parameters  were  obtained 
through  RRKM  calculations  (see  Appendix  and  Table  II)  using  a transition-state 
model  ( six-Piembered  ring)  . The  calculation  fits  the  data  for  the  following 
high-pressure  Arrhenius  parameters  (Figure  2) ; 

log  (kg/s-M  = 1200  - 38. 4/^  <T>  = 800  K (8) 

The  activation  energy  is  believed  to  fit  the  data  to  within  t 1 kcal/mol 
with  a concomitant  uncertainty  of  ±.  .27  in  log  (Ag/s“*). 

C . Benzylvinylether 

In  an  attempt  to  study  the  bond-breaking  process  to  produce  ethanal-2-yl 
radical  ("CHjCHO)  without  interference  from  a lower  energy  (molecular  elim- 
ination) reaction  path,  benzylvinylether  was  decomposed.  The  reaction  system 
under  investigation  was  the  following; 


C7H7OC2H3  — 

+ CH2CHO 

(9) 

k 

C^H^'  + DI  - C^H^D 

+ I’ 

(10) 

The  extent  of  decomposition  of  benzylvinylether  was  followed  by  trapping 
one  of  the  reaction  products,  benzyl  radical,  with  an  excess  of  DI  to  produce 
toluene-dj  (C^H^D) , which  could  readily  be  monitored  at  m/e  = 93.  This  method 
was  necessary  because  benzylvinylether  had  only  a very  small  intensity  molecular 
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Table  II 


MOLECULAR  PARAMETERS  ENTERING  THE  RRKM-CALCULATION  OK 
THE  RETRO-ENE  REACTION  OF  t-BUTYLVINYLETHER . t-C^HgOC^H,  - i-C^Hg  + CH3CHO 


1 

Molecule 

0 CH3 

11’ 

Activated 

Complex 

CH3 

H 

2940(6) 

2940(6) 

2900(3) 

2900(3) 

287013) 

2870(3) 

Frequencies 

1640(1) 

1640(1) 

1430(3) 

1430(3) 

and 

1380(5) 

1380(5) 

1220(4) 

1220(3) 

degeneracies 

1000(4) 

1000(4) 

900(3) 

900(3) 

650(4) 

650(4) 

580(1) 

580(1) 

1 

j 

450(4) 

450(4) 

1 

1 

350(4) 

350(4) 

• 

250(1) 

330(3) 

j 

220(3) 

250(1) 

1 

80(1) 

220(2) 

1 

i 

60(1) 

3 

10»2Oi  [ I / (gr  cm*) 

ABC 

2.05  * 10* 

2.05  ^ 10* 

E®  / kcal/mol 

38  .0 

' / e.u. 

1 

93  .40 

85  .80 

1 loglAgoo  / s-‘) 

1 

12 .00  1 

1 

1 

1 £^00  / kcal/mol 

38.4 

/ s-‘) 

11  .83 

1 

e”oo  / kcal/mol 

38  .0 

i 

•i 
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Ion  peak  m/e  = 134  (<  1%  of  base  peak  at  m/e  = 91)  and  other  minor  fragmen- 


! 

I 


tation  peaks,  so  that  it  proved  to  be  impossible  to  monitor  accurately  the 
disappearance  of  the  parent  compound  through  its  fragmentation  pattern.  The 
concentration  of  the  excess  component  D1  was  chosen  to  give  a yield  >95?! 
for  reaction  using  the  recently  measured  rate  constant  kj^  at  these  temperatures.® 
Special  attention  was  paid  to  the  requirement  of  working  at  low  concentrations 
of  benzylvinylether  ^ 3 * 10*“*  molec  s”*)  in  order  to  suppress  the 

recombination  reaction  of  benzyl  and  ethanal-2-yl  radicals  to  form  3-phenyl- 
propanol,  which  was  indeed  identified  by  its  mass  spectrum  (m/e  = 134,  105,  92, 

91)  at  flow  rates  of  benzylvinylether  higher  than  3 * 10*“*  molecules  s~*  . Figure  3 
displays  the  results  of  the  unimolecular  decomposition  of  ben  ylvinylether  at 
flow  rates  ^ 3 * 10*“*  molec  s~*  in  the  form  k^jj^j^/s"*  versus  T/K. 

In  order  to  determine  the  high-pressure  Arrhenius  parameters  for  the 
bond  scission  process,  RRKM  calculations  were  performed  using  the  well  known 
vibrational  model  often  mentioned  in  the  literature.®  (See  Appendix  and  Tables 
111  and  IV).  The  results  of  these  calculations  are  plotted  in  Figure  3,  together 
with  the  experimentally  determined  rate  constants  • The  calculation  fits 

the  data  for  the  following  high-pressure  Arrhenius  parameters: 

log(kg/s-‘)  = (15.05  ± .30)  - (49.30  i 1.00)/e  <T>  = 750  K (11) 
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THERMOCHEMiCAL  QUANTITIES  (Standard  State:  1 atm) 


Reference 


Table  IV 


MOLKCULAK  PAllAMKTKltS  KNTKKING  THK  KKKM-CAIX-'UUTION  OK 
Tia:  HATE  CONSTANT  t'OR  THE  KKIACTION  HENZYEVINYE  ETHEH  - C^H^  ’ + ‘011^0110 


Act Ival ed 

Molec ule 

^5v 

Cemple.\ 

s' 

» ' 

0 

L ' 

3080(7) 

3050(0) 

2030(3) 

2820(1) 

1500(5) 

1740(1) 

1300(5) 

1510(0) 

i 170(5) 

1330(3) 

1 100(1) 

1110(0) 

K requeue les 

1020(3) 

000(5) 

050(0) 

880(3) 

and 

830(3) 

700(1) 

000(2) 

000(2) 

degeneracies 

010(1) 

010(2) 

510(3) 

510(1) 

340(2) 

420(4) 

280(4) 

380(2) 

200(2) 

250(1) 

, 80(1) 

: 

150(4) 

100(1) 

r(C-O)/ A 

1 .33 

1.33  x 2.5 

3 .32 

a 

1 I / (gr 

ABC 

3 .055  lO** 

8 .050  » 10'* 

1 / I 

1 .45 

1 O'*  “ 1 ( pheny  1 ) / g r cm''* 

73  .0 

10'*“  1 ^ ( ben/.y  l-2-et  hana  lyl) 

50  .0 

1 ^JOO 

1 

00  .80 

105.02 

>“K  ''aoi.  ^ 

14  .50 

'‘•K  ''t*..  **■* 

15  .05 

K*'  kiM  1 tmtl 

17  .0  I 

t k.  t 1 Mill 

1 

1 

10  ..1 

1 11  (hi  I'liviii  I •*  % 


i 

[ 


f 

t 


I 


I 

) 


i 

f 

I 


f 
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D • Reactions  of  2-Ethanalyl  Radical  (‘CH..CHO) 

Two  groups  of  experiments  were  performed:  (1)  one  In  which  did 

not  exceed  approximately  8 * 10‘®  molecules  s"‘  ("low"  flow  rates  of  benzyl- 
vlnylether)  and  (2)  the  other  with  ^ 5 * 10**  molecules  s~‘  ("high" 

BVCi 

flow  rates)  . 

At  "low"  flow  rates  of  benzylvinylether,  CHjCHO  (m/e  = 44,  43,  29), 

CH^CO  (m/e  = 42,  14),  and  CO  (m/e  = 28)  were  detected  mass  spectrometrically . 

Upon  addition  of  DI , CHjDCHO  (m/e  = 45,  44,  29)  appeared,  CH^CO  disappeared,  and 
CO  formation  was  attenuated  with  respect  to  experiments  without  DI  . Therefore, 

It  is  believed  that  at  low  flow  rates,  the  reactions  (12)  and  (13)  are  operative: 


2CH2CHO 

CH3CHO  + CH2=C=0 

(12) 

CHgCHO 

- CH3CO  -*  CHj  + CO 

(13) 

The  fast  homogeneous  disproportionation  reaction  (12)  apparently  takes  place 

even  at  low  flow  rates,  whereas  the  decomposition  of  CHjjCHO  (13)  Is  thought 

to  take  place  on  the  walls  of  the  reaction  vessel,  the  first  step  being 

* 

exothermic  by  about  10  kcal/mol . We  attempted  direct  detection  of  the 
radicals  by  using  low  electron  voltage  (LEV) , but  the  Intensities  were 
below  our  limits  of  detectability. 

At  high  flow  rates  of  benzylvinylether,  LEV  experiments  indicated  no 

CHjCHO  radicals  present;  a finding  which  was  further  supported  by  the  results 

AsNUSitng  a stabi  1 1 'at  lim  energy  for  2-ethanalyl  radical  of  2 .(M)  kra  I/mol. 

See  tllaruaolon 
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of  oxperimentM  of  ilecomposl  t ion  of  bun '.ylvinyl  ether  with  DI  , which  gave 
tndlcutlon  of  the  presence  of  only  simxll  amounts  of  CH^UCHO  or  CH^IX^O 
(m  e If*).  However,  In  the  runs  with  1)1,  CH3CHO  and  CH^CO  were  found 
which  means  that  at  high  levels  of  CH^t'HO  the  disproportionation  reaction 
(12)  Is  suftlctenlly  fast  to  compete  with  reaction  (1-1): 

CH^flK)  + Dl  - CH^IX'HO  + T (It) 

Kurthenm*re,  In  rvms  with  1)1,  t’Hj  radical  was  Indirectly  detected  as 
t’HjD  (m  e 17,  It*,  IS)  and  t'O  occurrence  was  again  attenuated  due  to 
removal  »*f  CH^l’HO  radical.  A series  of  fast  secondary  reactions,  probably 
Involving  H-atoms,  apparently  takes  place  at  higher  flow  rates  as  Indicated 
by  the  appreciable  Intensities  of  the  peaks  at  m/e  = 78,  79,  lOt , 105,  106. 

IV  DISCUSSION 

The  high-pressure  Arrhenius  activation  enei*gy  for  the  ret  ro-ene  reaction 
of  ethy Iv Inyl ether  Is  in  good  agreement  with  the  values  reported  In  the 
literature.*  RRK  calculations**  suggest  that  the  occurrence  of  the  highei' 
enejgy  reaction  path  of  C-0  bond  fission  at  1300  K may  have  escaped  the 
present  detection  techniques  due  to  the  unfavorable  ratio  of  the  rate  constants 


Uiglk,  s-')  ^ 11.5  - 11.5/P,  loglk.^  s-')  15.0  - (65.0  t 2)  P,  where  k,  refers  j 

to  the  r»*t  ro-ene  react  u>n  anti  k^  refers  ti*  the  l‘-<)  lM>nd  sctsslt'n  rt*actton  | 
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of  the  molecular  elimination  mode  (kj)  to  the  homolytlc  C-O  bond-scl sslon 
mode  (kj).  This  argument  Is  further  strengthened  If  one  keeps  in  mind  that 
the  C-O  bond  strengths  are  known  to  be  accurate  only  to  + 2 kcal/mol.  In 
the  case  of  the  thermal  decomposition  of  t-butylvlnyl ether , RRK  calculations 
indicate  a ratio  ki/k2  = 220  at  T = 900  K,  which  almost  certainly  rendered 
the  detection  of  the  higher  energy  reaction  path  with  respect  to  the  retro-ene 
reaction  impossible  using  the  present  experimental  technique.  The  Arrhenius 
activation  energy  for  the  retro-ene  reaction  of  t-butylvlnyl ether  falls  within 
the  expected  range.  The  present  value  is  higher  and  believed  to  be  more 
reliable  than  the  one  from  the  classical  high-pressure  static  experiment** 

I for  reasons  cited  above  (narrow  range  of  temperatures) . Since  the  attempt 

i 

I of  studying  the  reactivity  of  ethanal-2-yl  radical  failed  using  ethylvinyl  ether 

and  t-butylvinyl  ether  as  a source,  we  chose  a compound  whose  lowest  energy 
reaction  path  was  a bond  fission  process  to  give  the  desired  radical  (CH^CHO)  . 
The  bond  fission  of  the  C-O  bond  in  benzylvinyl ether  is  the  lowest  energy 
reaction  path  and  yields  benzyl  and  ethanal-2-yl  radicals  unencumbered  by  the 
reaction  products  of  the  retro-ene  reaction.  The  determination  of  the 

Arrhenius  activation  energy  for  the  bond  split  reaction  in  BVE  yields  the  heat 

i 

of  formation  and  the  stabilization  energy  of  ethanal-2-yl  radical  under  the 

i 

^ log(k,  s"*)  = 11.85  - 38.0/0,  loglk^/s-*)  = 15.0  - 60  9 . 
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constraint  of  the  following  assumptions;  (1)  known  standard  heats  of  formation 
of  benzyl  radical  (AHf  (C_H-')  ),  (2)  known  C_,  values  for  benzylvinyl  ether , 

'^•300  r' 

benzyl,  and  ethanal-2-yl  radical  as  a function  of  temperature,  (3)  known 
overall  entropy  change  for  the  bond  fission  reaction  (AS°),  and  (4)  known 
rate  of  recombination  of  the  radicals  at  the  temperature  of  decomposition  of 
benzylvinyl ether . 


The  Arrhenius  activation  energy  for  the  homolytic  C-0  bond  fission  of 
benzylvinylether  may  be  identified  with  AE°  for  the  reaction  if  the  activ- 
ation energy  for  the  recombination  reaction  of  benzyl  and  ethanal-2-yl 
radical  is  assumed  to  be  zero  at  the  temperature  of  decomposition.  Accordingly, 
AH°  = 50 .8  ±.  1 .0  kcal/mol  at  750  K and  with  the  C values  of  Table  III, 

P ,500 

one  obtains  the  bond  dissociation  energy  ( = AKg^^)  at  300  K for  reaction  (9) : 
(9):  ^5oo  = 52 .2  ± 1 .0  kcal/mol.  The  heat  of  formation  of  the  species 
involved  in  reaction  (9)  are  related  through  relation  (15)  : 


BDECC^H^-OCgHg) 


= ^300  = AH° 


>300 


,300 


(CH^CHO)  - 


- AH°  (C.H.OC.H,) 

f,300  ^ ^ 2 ^ 


(15) 


With  the  group  additivity  value  of  —1.73  kcal/mol®  for  the  standard  heat  of 
formation  of  benzylvinylether,  and  a recent  value  for  the  heat  of  formation 
fo  benzyl  radical  ( 46.5  kcal/mol*  ),  one  obtains  a value  of  3.90  +2.0 
kcal/mol  for  (CH„CH0)  . On  the  basis  of  a value  of  98.0  kcal/mol  for 
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the  bond  strength  in  H-CH^CHO , corresponding  to  BDE  of  a bond  in  a nornuil 
alkane,  (CH„CHO)  would  be  7.1  kcal/mol . The  difference  between  the 

"hypothetijpal"  standard  heat  of  formation  of  CH^CHO  (7.1  kcal/mol)  and  the 
above  experimental  value  of  3.9  ±2.0  kcal/mol  is  the  resonance  stabilization 
energy  (SE)  due  to  delocalization  of  the  unpaired  electron  over  the  3-center 
Tf-system  analogous  to  allyl  radical  and  is  computed  to  be  3 .2  ± 2 .0  kcal/mol. 

This  value  for  the  resonance  stabilization  energy  compares  well  (within  the 
uncertainties)  with  the  values  recently  found  for  the  closely  related  acetonyl 
radical  (CBjCOCH^)**  and  is  significantly  lower  than  the  allyl  resonance 
stabilization  energy  of  12.0+1  kcal/mol.®  If  the  activation  enei-gy  for  the 
C-0  bond  fission  of  benzylvinylether  is  identified  with  AH®  for  the  reaction 
(i.e.,  the  activation  energy  for  the  recombination  of  benzyl  and  ethanal-2-yl 
radical  is  assumed  to  be  zero  at  0 K , and  RT  at  the  temperature  of  decomposition), 
one  computes  AH°  (CH„CH0)  = 2.4  ± 2.0  kcal/mol  and  SE  = 4 .7  + 2 .0  kcal/mol. 

* »300  ^ 
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Figure  1 Unimoiecular  Decomposition  Rates  as  a Function  of 


Temperature  for  the  Molecular  Elimination  Mode  of  Ethylvinylether . 

( A . • > O represent  different  series  of  experiments) 

Figure  2 Unimoiecular  Decomposition  Rates  (k  ./s"‘)  as  a Function  of 

. . urxi. 

Temperature  for  the  Molecular  Elimination  Mode  of  . -Butylvinyether . 
The  loss  of  t-Butylvinylether  was  monitored  at  m/e  = 57  in  series 
o,  and  at  m/e  = 100  in  sex'ies  •)  . 

Figure  3 Unimoiecular  Decomposition  Rates  (k  ./s“‘)  as  a Function  of 

llUX 

Temperature  for  the  Bond  Fission  in  Benzylvinylether . 

( •,  high  flow  rate  series;  o,  low  flow  rate  series;  see  text 
for  details)  . The  dashed  line  corresponds  to  E ±1  kcal/mol 
at  <T  > = 750  K. 


Appendix 


f 

I 

I • In  this  section,  details  about  the  RRKM  calculations  shall  be  given. 

i 

! The  reader  is  referred  to  Tables  I-IV  of  the  text  which  display  the  actual 

numerical  values  used  in  the  calculations. 


Ethylvinylether 


A rather  complete  structural  investigation  of  methylvinylether  (C2H3OCH3)* 
indicated  the  existence  of  nonbonded  attractive  interactions  between  the 
methyl  and  vinyl  groups,  a result  recently  obtained  by  ab  initio  calculations 
as  well.*  Calculating  the  entropy  of  methylvinylether  at  300  K,  following 
group  additivity  methods,*  yields  (CH3OC3H3 ) = 73.50  e.u.,  whereas  the 

r;.'sult  of  an  entropy  calculation  on  the  basis  of  the  assigned  frequencies 
experimental  geometry  (cis-planar  conformation)  of  methylvinylether*  reveals 
a value  of  69.30  e.u.,  reflecting  the  nonbonded  attractive  interactions  or 
higher  degree  of  stiffness  in  that  molecule.  About  3.0  e.u.  of  the  total 
entropy  difference  of  4.20  e.u.  can  be  accounted  for  by  assigning  substantial 
barriers  to  internal  rotation  to  the  0-vinyl  intenial  rotation  (3.3  kcal/mol) 
and  0-methyl  internal  rotation  (3.0  kcal/mol)  in  methylvinylether.  The  group 
additivity  value  for  the  entropy  of  ethylvinylether  is  83.30  e.u.  and  was 
corrected  by  3.0  entropy  units  to  give  80.30  e.u.  to  account  for  the  barriers 
to  internal  rotation.  Except  for  the  Introduction  of  CHj * Instead  of  H’ , only 


i 


*A  barrier  to  internal  rotation  of  about  3.3  kcal/mol  for  the  0-vinyl  torsion 
results  in  the  same  entropy  contribution  as  the  experimentally  observed  torsion 
of  frequency  v = 220  + 20  cm“*  (reference  l) . 
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minor  corrections  of  the  vibrational  frequencies  of  methylvinylether  had 
to  be  carried  out  in  order  to  predict  the  frequencies  for  ethy Ivinylether 
to  give  a value  of  (C^HgOCgHg ) = 80.30  e.u.  The  frequencies  of  the 
six-membered  transition  state  for  reaction  (2)  were  then  chosen  in  such 
a way  to  result  in  the  experimentally  observed  Arrhenius  A-factor  of  lO**’^’' 
s"*  at  <T>  = 900  K (the  temperature  dependence  of  log  A is  negligible, 
c.f.,  Table  I).  Basically,  the  three  hindered  internal  rotations,  together 
with  some  of  the  low-frequency  bending  modes  of  ethy Ivinylether  were 
replaced  by  higher  frequency  out -of -plane  ring  bending  modes,  which  were 
only  slightly  lower  than  the  corresponding  out-of-plane  frequencies  in 
benzene,  in  order  to  simulate  the  stiff  six-membered  ring  transition  state, 
where  six  electrons  were  redistributed  similar  to  an  aromatic  hydrocarbon 
^-system.  A list  of  the  molecular  parameters  for  the  molecule  and 
activation  complex  is  presented  in  Table  I of  the  text. 

t-Butylvinylether 

The  published  high-pressure  Arrhenius  A-factor  for  reaction  (6)^  was 
recognized  to  be  too  low  together  with  too  small  a value  for  the  activation 
energy,  probably  due  to  the  rather  restricted  range  of  temperatures  used  in 
assessing  the  Arrhenius  parameters.  The  A-factor  was,  therefore,  obtained 
by  making  a correction  of  + .48  to  log  A = 11 .47  for  the  corresponding- 
reaction  of  ethylvinylether  due  to  reaction  path  degeneracy  of  3,  thereby 
obtaining  log  A -S'  12.00  for  reaction  (6)  at  <T>  = 800  K.  The  concomitant 
assumption  is  that  with  respect  to  ethylvinylether,  the  lower  frequency 
torsions  due  to  the  heavier  masses  in  t-butylvinylether  are  balanced  by 
lower  ring  puckering  or  ring  out-of-plane  deformation  modes  in  the  .activated 
complex.^  The  group  additivity  value  for  t-butylvinylether  (SJqo  ~ 96.40  e.u.) 
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was  decreased  by  3.0  e.u.  to  give  93.40  e.u.  to  account  for  (assumed)  attractive 
nonbonded  interactions  analogous  to  the  case  of  ethylvinylether . The  frequen- 
cies of  the  molecule  were  chosen  to  correspond  to  S3QQ  = 93.40  e.u.,  and  the 
ones  for  the  activated  complex  were  determined  in  such  a way  to  be  compatible 
with  log  A = 12.00.  Table  II  discloses  the  details  of  the  molecule  and 
activated  complex  parameters  entering  the  RRKM  calculation. 


Benzylvinylether 

The  Arrhenius  A-factor  for  reaction  (9)  was  computed  from  the  overall 
entropy  change,  AS°  (Table  III),  and  the  assumed  rate  constant  for  the 
recombination  reaction  (-9),  logdc.g/M"*  s"^)  = 10.00.  The  assumption 
concerning  the  recombination  reaction  in  the  present  study  is,  that  it  has 
no  activation  energy  at  T = 750  K (Ej.  = 0) , so  that  one  can  identify  the 

Arrhenius  activation  energy  E^^  for  the  bond  fission  process  (9)  with  AE° 

for  the  reaction.  This  yields  the  following  A^j-factor  at  T = 750  K through 
relation  (16) : 

10 1 o * 

Rln  — = -33.12  +10,16  log(A^/s-‘)  = 15.02  (16) 

The  value  of  15.02  for  log  A^j  compares  well  with  the  ones  for  the  similar 
bond  scission  processes  (17)  and  (18) : 

CgHsC^Hg  + CH3  log(Ajj/s-‘)  = 14.60  at  <T>  = 950  K ^ (17) 

CeH5C3H7  - C7H7'  + CgHg  ‘ log(Ajj/s“‘)  = 14.90  at  <T>  = 950  K ® (18) 


* / 

10,16  represents  the  correction  term  (R(l  + InR  T^)  ) for  the  change  of 

standard  states;  R is  the  gas  constant  in  units  of  cal  mol“‘  1 and 

R^  is  in  units  of  -t  atm  mol"^  K~‘  . 
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The  group  additivity  value  for  Sj,,^  of  benzylvlnylether  was  taken  to 
be  too  large  by  ~3  cal  K“‘  mol“‘  as  in  the  case  of  methylvlny lether  for 
reasons  of  intramolecular  interactions,*  so  that  the  molecular  vibrational 
frequencies  and  moments  of  Inertia  were  chosen  to  match  = 99.80,  instead 

of  SjQQ  = 103.60.®  The  entropy  of  the  activated  complex  was  then  chosen  to 

be  compatible  with  log  = 15.02.  Table  IV  shows  the  details  of  the  I 

molecular  and  activated  complex  parameters  for  the  RKKM  calculation. 
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Chapter  4 


Infrared  Photodecomposition  of  Kthyl  Vinyl  Ether. 

A Chemical  Probe  of  Multiphoton  Dynamics 

Sir; 

Chcmi.siry  following  niullipholon  ab.sorption'  may,  in 
general,  involve  eollisionless  and  collision-induced  processes. - 
The  d\  namics  and  time  regimes  associated  with  these  proces.ses 
remain  to  be  well  characterized  experimentally.  We  report  here 
the  COj  TEA-lascr  photolysis  of  ethyl  vinyl  ether  (EVE.  1) 
and  its  relevance  to  these  questions. 

Conventional  pyrolysis'  of  EVE  yields  ethylene  and  acet- 
aldehyde by  a retro-enc  molecular  elimination  (eq  1 ). 


Ht  CH  . 

CH, 

■i  1 

1 

. ir 

1 C CH, 

CH, 

CH, 

/' 

it 

log.k,  s”* 

= n..i  - (1, 

'Z.203R'T 

Decomposition  with  a focu.sed  laser  yields  not  only  the.se 
conventional  products  but  also  comparable  amounts  of  kctenc, 
ethane,  and  butane.'*  This  suggests  that  a simple  bond  fission 
process'  (cq  2)  competes  effectively  with  eq  I upon  multi- 

h'  ^ch. 
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photon  excitation.*’  Kctene  is  formed  by  subsequent  dispro- 
portionation reactions  of  -CHyCHO  and  the  ratio 
|CT1  iCH0|/[CH  'C01  provides  a measure  of  the  relative  rates 
of  processes  I and  2. 1 his  ratio,  as  a function  of  pressure,  is  as 
follows:  1.89 (440 Torr).  1.89(280),  l.64(2.J0),  1.45(25).  1.44 
(20),  1.71  (1 1),  1.79(10).  1.86(10).  1.84  (5),  and  1.75  (10 -h 
40  torr  of  He).  Thus,  k \ ^ 0.37  over  the  pressure  range 
studied.  Irradiation  with  an  imfocu.sed  beam  yields  only  eth- 
ylene and  acetaldehyde. 

The  two  reactions  compete  at  high  energies  since  the  reac- 
tion channel  density  for  eq  2 exceeds  that  for  eq  I (as  reflected 
in  their  .4  factors).  A number  of  molecules  will  exhibit  similar 
behavior  if  their  lowest  thermal  path  involves  a cyclic  transition 


state,*’  but  FvVE  is  especially  interesting  since  the  difference 
in  activation  energies  for  its  two  lowest  energy  channels  is  large, 
thereby  increasing  the  dynamic  range  available  for  studying 
the  energy  distribution  of  reacting  molecules. 

Our  experimental  results  can  be  considered  in  terms  of  the 
following  steps:  laser  pumping.  A(E)  -*  A(E')  |Ap(l,E)l;col- 
lisional  energy  pooling  or  deactivation.  M + A(E')  *54  + 
A(E''),  (AlI;  and,  reaction  via  channel  i,  A(E)  -♦  (products),, 
lilt,(E)j.  Such  a scheme  calls  attention  to  the  various  competitive 
processes,  each  of  which  may  dominate  under  different  ex- 
perimental conditions.’ 

In  our  pressure  range.  ~1  100  collisions  occur  during  the 
la.ser  pulse.  The  observed  lack  of  pressure  dependence  thus 
suggests  that  most  the  chemistry  occurs  after  the  pulse  Fur- 
thermore. assuming  the  applicability’-'*"’  of  a quantum  RRK 
model,"  we  find  that  A |(E)  and  A- 2(E)  are  comparable  at 
~10*’  10’  .s~',  which  is  < collision  frequency,  imply  ing  that 
our  results  are  predominantly  collisional.  Collisions  between 
molecules  in  the  irradiated  region  occur  rapidly  compared  with 
c.scape  and  sub.sequent  cooling,  since  the  mean  free  path  is 
small  at  our  pressures.  (Isotopic  specificity  may  still  obtain 
under  these  conditions  since  collisions  with  "cold"  molecules 
arc  simply  deactivating.)  Therefore,  in  our  experiments,  the 
chemistry  appears  to  be  that  of  a collisionally  (V  V)  thcr- 
malizcd  system.  This  is  also  consistent  with  recent  experi- 
ments'*’which  suggest  that,  at  energy  flucnces  of  ~1  Jcm~’. 
although  considerable  excitation  occurs,  only  a small  fraction 
of  excited  molecules  react.  The  low  conversions  ty  pically  ob- 
served in  this  regime  are  thus  consistent  with  an  energy  dis- 
tribution function  that  decreases  monotonically  above  £.\..  the 
threshold  for  reaction.  Consequently,  the  collisional  nature  of 
this  chemistry  indicates  the  Arrhenius  forms  of  k 1 and  k : can 
be  used  in  expressing  the  yield  ratio,  giving  a temper, tture  of 
~1600K.Thc  tempcraturedcpcndcnceof  A'l/A:  (for  7 = 300. 
1000.  2000  K,  (( i//t;  = (0.55  X 10'’,  18.94,  8.68  X lO"’)  in- 
dicates the  sensitivity  of  the  branching  ratio  of  EVE  as  a probe 
of  excitation.  The  cleanliness  of  IR  laser  photolysis  relative 
to  the  high  temperature  pyrolysis’  of  EVE  suggests  this  is  an 
attractive  alternative  to  conventional  activation  techniques. 
Evidence  indicates  that  increasing  fluence  will  raise  the  ap- 
parent “temperature"  ' ' so  that  eq  2 would  compete  more  ef- 
fectively with  cq  1. 

In  summary,  we  have  developed  a sensitive  chemical  probe 
of  energy  distribution  applicable  to  bulk  and  collisionless 
systems.  Our  results  are  consistent  with  recent  studies'"  in 
terms  of  lifetimes  and  energy  distribution  functions.  We  are 
currently  studying  the  intensity  and  buffer  dependence  of  k i/A  2 


1 


i 

i 

4 


Ji 


and  extending  our  experiments  to  include  the  collisionless  re- 
gime by  increasing  fluence  and  decreasing  pressure. 
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Chapter  5 


APPLICATION  OF  RRKM  THEORY  TO  THE  REACTICWS 
OH  + NO^  - HONOj,  + Nj,  (1) 

and 

CIO  + NOjj  + Njj  - ClONOj  + Njj  (2); 

A MODIFIED  GORIN  MODEL  TRANSITION  STATE. 

G.  P.  Smith  and  D.  M.  Golden 
Thermochemistry  and  Chemical  Kinetics  Group 
Stanford  Research  Institute  , Menlo  Park,  California  94025 

ABSTRACT 

We  have  calculated  rate  constants  as  a function  of  both  temperature 
and  pressure  for  the  title  reactions  using  RRKM  theory  In  conjunction  with 
a modified  Gorin  transition  state.  The  modification  Introduces  a hindrance 
parameter  which  accounts  for  repulsive  Interactions  between  the  rotating 
fragments . 

At  the  highest  stratospheric  pressures  (~  50  torr)  and  at  stratos- 
pheric temperature  ( ~ 220°K)  , the  extent  of  "fall-off"  from  first-order 
[Nj]  dependence  Is  ~ 70%  for  reaction  (1)  and  ~ 35%  for  reaction  (2) . 
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INTRODUCTION 


Nitric  acid  formation  via, 

OH  + NO2  + Nj  - HONO^  + Ng  (1) 

is  an  important  stratospheric  reaction  which  functions  as  a sink  for  ozone 
destroying  HO^  and  NO  cycles.*  Thus  rate  constants  are  available  at 

^ X 

several  temperatures  and  nitrogen  pressures.*"®  At  stratospheric  pressures, 
reaction  (1)  is  in  the  "fall-off"  region,  between  the  low-pressure  termolec- 
ular  limit  and  the  high-pressure  limit  where  bimolecular  formation  of 
metastable  nitric  acid  is  rate  determining.  Hence,  reaction  (1)  is  an 
excellent  candidate  for  RRKM  calculations®  in  order  to:  (a)  see  if  the 
data  can  be  fit  within  the  confines  of  the  theory,  using  reasonable  model 
parameters;  (b)  provide  a basis  for  extrapolation  and  interpolation  of 
reliable  reaction  (1)  data  to  other  temperatures  and  pressures;  and  (c) 
pursue  a priori  criteria  for  the  future  selection  of  transition  state 
parameters  in  simple  bond  fission  reactions. 

In  furtherance  of  this  last  goal,  we  have  adapted  our  successful  model 
for  reaction  (1)  to  predict  stratospheric  rates  of  the  analogous  reaction: 

CIO  + NOjj  + Nj  - ClONOj  + Njj  (2) 

This  reaction  is  being  extensively  examined  as  a sink  for  ozone-destroying 
stratospheric  chlorine  species,  which  are  introduced  via  chlorofluorocarbon 
propelled  aerosols,'^  but  very  few  experimental  rate  constants  are  currently 
available . ® 
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Recent  advances  in  the  understanding  of  unimolecular  reactions  have 
made  it  clear  that  a simple  "fixed"  transition  state  is  an  inadequate  model 
for  simple  bond  fission  reactions.  These  advances  have  come  from  two 
directions,  the  microcanonical  approach  as  exemplified  by  the  work  of  Troe 
and  coworkers,®  and  the  canonical  approach  which  we  have  favored.^® 

In  the  microcanonical  view.  Quack  and  Tro^  have  illustrated  the 
problem  by  pointing  out  that  individual  open  reaction  channels  may  have 
maxima  at  various  locations  along  a specified  reaction  coordinate.  They 
have  formulated  their  "Adiabatic  Channel  Model,"  taking  this  into  account 
and  using  a single  universal  parameter  to  describe  the  potential  surfaces 
for  bond  scission  in  a number  of  small  polyatomic  systems.  (They  have  also 
recently  presented^c  a canonical  model  which  in  some  ways  is  similar  to  that 
to  be  presented  herein.) 

In  the  canonical  view,'Owe  have  postulated  a transition  state  which 
can  become  "tighter"  as  the  temperature  increases.  As  a result  of  the 
decrease  in  position  of  the  centrifugal  maximum,  the  rotation  of  the  product 
fragments  is  increasingly  restricted.  We  thus  introduce  a hindrance 
parameter,  Tl,  to  describe  this  tightness.  If  Tj  is  independent  of  temperature. 


we  have  a model  which  is  equivalent  to  Reference  9c;  if  it  is  T dependent,  we 
have  introduced  a second  parameter. 

The  recombination  rate  kj  can  be  deduced  from  the  RRKM  unimolecular 
decomposition  rate  k-j , since  kj  = Kgqk_j . RRKM  theory®  is  a statistical 
theory  which  assumes  fast  random  access  by  allowed  Interaal  states  to  molecular 
energy  E* , and  a unique  critical  configuration  along  one  coordinate,  which  is 


intermediate  to  all  decompositions.  Hence  the  Important  parameters  for 
the  calculation  relate  to  the  structure  of  this  transition  state.  The 
rate  is  given  by ; 

Q+  e-Eo/kT  r G(E+)e-E''/kT  dE+ 

k . = / (3) 

hQjQj  / 1 QtG(E+)/hQ,N(E*)F3a) 

where  E^,  is  the  critical  energy,  E is  the  total  energy,  and  E = E ~ E^ 
is  the  maximum  energy  of  the  critical  configuration;  G(E"'’)  is  the  sum  of 
the  states  of  the  complex  below  energy  E"*"  and  N(E*)  is  the  density  of 
molecular  states;  Q^/Qi  is  the  partition  function  ratio  of  the  inactive 
modes  (moments  of  adiabatic  rotors) , and  is  the  partition  function  for 
the  active  molecular  modes;  F is  the  Waage-Rabinovitch^ * centrifugal  correc- 
tion term  for  conservation  of  angular  momentum  among  the  adiabatic  rotations, 
(JO  is  the  (Lennard-Jones)  collision  frequency  of  excited  substrate  with  bath 
gas  at  pressure  P,  and  3 is  the  collisional  efficiency  for  stabilization 
(sufficient  energy  removal).  Our  calculations  used  the  Stein-Rabinovitch 
algorithm**  to  compute  G(E'*’)  and  N(E*)  . 


DETAILS  OF  THE  CALCULATION 

For  our  calculations,  is  computed  from  the  JANAF  Tables,*®  These 
same  tables  furnished  the  molecular  frequencies  for  computing  N(E  ),  and 
the  dissociation  energy  of  HONO^  for  computing  Eg.  We  chose  the  top  of 
the  centrifugal  barrier  for  a Lennard-Jones  potential  as  the  position  of 
the  critical  configuration  along  the  bond  axis;  we  have*'* 

Q^/Q,  = I'^/l  = (r'*’/ro)*  = (SD^/RT)*'^^  (4) 

where  r is  the  distance  between  the  OH  and  NO^  centers  of  mass  and  and 
I are  moments  of  inertia  in  the  transition  state  and  molecule  for  the  two- 
dimensional  external  rotation  with  axis  perpendicular  to  the  HO-NO^  bond. 
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The  collision  frequency,  u),  can  be  calculated  by  averaging  the  Lennard-Jones 

collision  diameters,  5.2  X for  HONO^  (estimated  using  the  viscosity-derived 

hard  sphere  diameter**  for  NO^Cl,  since  no  critical  data  exists  for  HONO^) 

0 

with  the  value  of  3.8  A for  Njj.  Only  the  specification  of  frequencies  in 
the  critical  configuration,  to  calculate  GiE"*")  , remains. 

Previous  RRKM  calculations*®  for  reaction  (1)  estimated  the  nitric  acid 
transition  state  vibrational  frequencies.  This  approach,  however,  is  further 
complicated  by  the  temperature  dependence  of  these  frequencies,  since  the 
critical  bond  stretching  distance  varies  (see  equation  (4)).  We  have  chosen  a 
modified  Qorin  model®  >*'*  for  the  transition  state.  In  the  Gorin  model,  the 
internal  modes  of  the  transition  state  are  simply  the  vibrations  and  rotations 
of  the  independent  OH  and  NO^  fragments  (again,  from  the  JANAF  Tables**)  . Thus 
we  have  here  five  internal  rotations,  corresponding  to  the  external  rotations 
of  the  fragments.  This  effectively  activates  the  external  HONOj  rotation  about 

the  HO NOg  bond  axis,  which  is  mostly  an  NO2  rotation,  and  leaves  only 

the  other  two  external  rotations  inactive  and  adiabatic  in  the  molecule  (i.e., 
these  modes  do  not  share  in  the  randomization  of  energy).  Since  the  internal, 
fragment  rotations  are  not  actually  free  at  the  OH-NO^  distance  of  the  critical 
configuration,  our  modification  introduces  the  aforementioned  hindrance  parameter, 
which  decreases  the  entropy  of  internal  rotation  by  decreasing  the  effective 
moment  of  inertia  of  the  rotor.  We  thus  decrease  the  number  of  available 
rotational  states  in  accordance  with  the  volume  of  rotational  phase  space 
excluded  because  of  the  other  fragment's  presence.  Such  a rotational  model, 
as  opposed  to  the  usual  vibrational  ones,  should  be  reasonable  for  recombination 
reactions  with  only  centrifugal  barriers,  because  they  have  loose,  distant 


transition  states. 


In  thest'  calculations,  the  tMO-ilimensiona  1 OH  rotor  amt  a two-dimensional 
SO.  external  rot  a t ton, chosen  to  exclude  rotation  about  the  bond  axis,  are  eacli 

hindered.  Thus,  Roll  , wliere  »‘ach  1 * - , 1—  ^ . Anv  model  for  estlmat- 

OH  NO,  • 

Inn  the  hindrance  of  the  OH  rotor,  fv'r  example,  actually  estimates  a value  for 
,1  — * ' . The  hindrance  does  not  apply  to  the  remaining,  one-dlmensional  NO, 
rotor,  which  corresponds  to  ttie  active  molecular  external  rotation.  Tins 
model,  as  av'plled  to  nitric  acid.  Is  Illustrated  In  the  appendix. 

We  thus  have  a model  with  two  unknv'wn  parameters,  the  hindrance  T]  and  tlu> 
collision  efficiency  There  are  some  reasonable  limiting  criteria  for  the 
values  and  temperature  dependences  of  these  variables  wlitch  are  acceptable  in 

f 

attempting;  to  tit  the  data.  Previous  experimental  low-pressure  rale  constants 
Indicate  the  relative  efficiencies  of  various  third  bodies  for  a given  recom- 
bination reaction.'"  For  N.,  we  expect  ^ ~ .3  to  .5.  Furthermore,  there 

is  some  low  pressure,  JHH'’  data'*  which  constrains  the  choice  of  6,  us  expected, 
to  values  neai'  .5.  Troe*'  has  deieloped  a theory  relating  and  I?  to 
the  average  energy  transferred  vier  collision.  Since  <,\1>  Is  not  expected  to 
change  much  with  temperature,  the  values  of  derived  from  I?  used  for  tlie 

fitting  should  show  only  limited  variation. 

It  Is  more  difficult  to  set  reasonable  limits  on  T.  By  taking  van  der 
Waals  radii  for  atomic  sizes  and  constructing  a miithemat  leal  "brick  wall" 
half  way  between  the  transition  state  fragments,  one  value  for  the  hindrance 
can  be  estinwted  by  the  solid  angle  excluded  to  each  rotor  by  the  wall.  Sucti 
a nKidel  gives  ~ -IS*”!-  for  2!UP' , but  fails  to  Include  the  effect  of  correlated 
motions  among  rotors.  Furthermore,  lai'ger  hindrances  are  required  before  one 
approaches  the  usual  re;:ime  of  loose  vibrational  transitu'n  states.  For  example, 
a I'n'i  hindrance  of  the  2 two-dimensional  Internal  rotors  ITB'l  each)  would 
correspond  in  entropy  to  four  BO-cm-’  vibrations  i300'’K)  . In  view  of  the 
difficulty  of  choosing  a priori  values  for  v,  for  frequencies  for  the  alternative 
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vibrational  transition  state  model),  the  goal  of  investigating,  thi 
of  T)  In  this  conceptually  simple  model  by  fitting  data  bet  om*-*  app^  ' 
Once  the  data  are  encoded  in  terms  of  this  simple  ",  parameter,  a »• .. 
transfer  to  similar  system^’  may  be  considered.  This  T)  param*-trt 
a simple  empirical  measure  of  "tightness"  within  the  conttM  »*t  a . 
type  model.  This  Gorin  model  has  the  advantage  that  the  lo»ei  ht-i 
of  the  transition  state,  as  compared  to  a vibrational  iiK>del,  i»  •• 
with  observation. 
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RESULTS 


Figures  1 and  2 show  the  best  RRKM  matches  to  the  data  at  296°K , and 
other  temperatures  between  220°  and  550°.  A good  fit,  within  experimental 
error,  is  obtained,  although  a 15%  variation  in  3 or  (1  - T])  generally 
still  fits  the  data.  This,  then,  is  roughly  the  uncertainty  in  our  RRKM 
extrapolation  of  the  falloff  region  data  to  the  high  pressure  (k„)  and 
low-pressure  (kj,)  rate  constants.  We  note  that  our  values  at  times  do  not 
agree  with  those  extrapolated  by  Anastasi  and  Smith^  using  Tree’s^®  Kassel 
integral  method  . 

Table  I lists  the  fitting  parameters  used  for  Anastasi  and  Smith's* 
data  at  220°-550°K,  and  Glanzer  and  Tree's  data*®  for  nitric  acid  decom- 
position at  1000°K.  The  1000°  "R  and  p are  consistent  with  extrapolation 
of  the  lower  temperature  trends.  The  fit  is  less  accurate,  but  the  data, 
and  hence  the  fitting  parameters,  are  more  uncertain. 

Table  I also  gives  values  for  <AE>  averaging  1.4  kcal/mole.  The  value 
of  <AE>  -^1.4  is  typical  of  nitrogen  energy  removal  rates  from  small 
molecules  such  as  CHjCN.*^  Although  considerable  latitude  exists  in  the 
choices  of  <^>,  the  table  indicates  <6E>  may  decline  slightly  with  tem- 
perature. (It  should  also  be  noted  that  values  of  3 are  strongly  dependent 
on  the  exact  calculation  of  the  "strong  collision  rate  constant."*’  Using 
Troe's  method  we  calculate  that  at  296°K,  3 ^ .37  as  opposed  to  .64,  as 
reported  in  Table  I.  This  difference  is  due  to  an  approximate  correction 
for  anharmonicity  and  a more  refined  treatment  of  the  torsional  motion,  and 
changes  <AE>  from  ~ 2 kcal/mole  to  ~ 0 .6  kcal/mole.) 

The  hindrance  reaches  unexpedtedly  high  values  at  moderate  temperatures, 
although  increasing  with  temperature  as  predicted  (since  r*  (eq.  4)  decreases 
with  temperature.  Internal  rotations  become  more  hindered) . This  suggests 
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that  the  netherland  between  hindered  rotations  and  torsional  vibrations  occurs 
at  distances  longer  than  our  original  model  envisioned.  (A  99*’  hindered 
rotation  mlKht  more  properly  be  viewed  us  a vibration.)  The  table  also 
indicates  some  temperature-to-temperature  irregularities  in  the  variation 
of  ^ and  T),  due  perhaps  to  ttie  sparsity  of  the  data  and  possible  errors  at 
some  t empe ra  t u res . 

The  2110"  entry  in  the  table  also  Illustrates  an  Identical  fit  usliii; 
different  parameters.  If  T)  is  lowered,  l^/I  must  be  treated  as  a new  variable 
and  decreased  to  maintain  the  same  entropy,  or  A-factor  k (I  - r,)(l^  1)]. 
However,  since  the  K-factor  (equation  3)  also  depends  on  I^/l,  P must  be 
Increased  to  balance  the  I^/l  decrease.  Such  adjustments  cannot  attain  more 
reasonable  hlKh-temperature  hindrance  values  and  still  retain  reasonable 
1^/1  and  [3  values. 

The  table  also  contains  two  alternate  fits  for  29(5",  which  depend  some- 
what upon  one's  interpretation  of  the  low-pressure  data.  Howard  and  Kvenson'* 
interpret  their  data  in  terms  of  an  underlying;  bimolecular  wall  reaction.  It  is 
sujjKested  from  our  calculations,  however,  that  their  experiments  are  not  in 
the  low-pressure  limit,  althoukh  both  experiments  and  calculation  are  nearly 
linear  in  this  region.  Thus,  great  care  must  be  used  in  interpreting  low 
pressure  results.  We  can  fit  the  data,  adjusted  for  their  wall  reaction 
(Figure  1,0=  .-18),  or  assuming  no  wall  reaction  (P  = .8-1),  altliough  the 

former  model  fits  the  data  of  refert'nce  3 better. 

Finally,  Figure  3 utilises  our  fit  of  the  data  and  the  standard  profile 
of  the  atmosphere'*"  to  calculate  (roughly)  the  eftective  bimolecular  rate 
constant  for  reaction  (1)  us  a function  of  altitude. 
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EXTRAPOLATION  OF  THE  MODEL  TO  CIONO^ 

Since  only  a few  rate  constant  measurements,®  all  below  10  torr, 
have  been  reported  for  reaction  (2) , the  necessary  inquiry  into  the 
degree  of  fall-off  from  third-order  kinetics  at  higher  stratospheric 
pressures  suggests  an  RRKM  calculation.  The  model  and  parameters 
chosen  should  be  consistent  with  the  successful  reaction  (1)  fitting. 

In  this  way  the  hindrance  concept  and  the  data  fitting  approach  permit 
an  educated  choice  for  the  reaction  (2)  transition  state. 

For  our  reaction  (2)  calculations,  (for  calculating  K„_) , the 
frequencies,  and  the  moments  of  inertia  for  the  CIONO^  molecule  were 
derived  from  Miller  et  al.,*'®  using  Dq  = 24.8  kcal/mole  The  corres- 

ponding JANAF  values^®  were  used  for  the  fragments.  The  critical  distance 
is  calculated  from  (4)  as  before. 

A Lennard-Jones  collision  diameter  of  ~ 5.6  X was  estimated  by 

adjusting  the  nitric  acid  value  for  the  larger  Cl  atom  size.  Since  true 

low-pressure  rate  constants  are  available  for  reaction  (2) , the  collislonal 

efficiency  is  a completely  empirical  value.  The  three  data  sets®  are  in  close 

s b 

agreement,  and  we  have  used  the  values  of  Zahnlser,  Chang,  and  Kaufman  in 
determining  P.  The  "efficiency"  is  thus  .83  at  298°  and  1.4  at  220°. 

While  this  temperature  dependence  is  consistent  with  that  using  He  as  a 
third  body,®^  and  with  the  reaction  (1)  fit,  the  magnitude  is  surprisingly 
large.  Such  a large  increase  in  P or  the  collision  diameter  upon  slight 
alteration  of  the  excited  molecular  species  would  be  very  surprising. 
Application  of  Troe's*’’  method  to  ClONOg-Nj  collisions  at  298°K  indicates 
that  P might  be  lowered  to  0 .69  as  a result  of  anharmonlcity  corrections 
and  the  effect  of  treating  the  internal  rotation  as  a torsion.  But,  this 
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would  still  mean  that  <AE>  ~ 2.8  kcal/mole  (compared  to  .6  kcal/mole  using 
Troe's  method  for  HONO^).  The  220°K  value  of  0 = 1.39  simply  reflects  some 
of  the  general  uncertainty  in  computing  the  strong  collision  rate  constant. 

It  should  also  be  kept  in  mind  that  if  E^,  were  slightly  higher  (1-2  kcal/mole), 
the  difference  in  0 and  <^>  between  HNO^-Nj  and  CIONO^-Nj  would  disappear. 

The  hindrance  parameter  was  estimated  in  the  following  manner.  We 
applied  our  simple  "brick  wall"  van  der  Waals  model  to  the  CIONO^  transition 
state  at  r'*’ , to  obtain  a predicted  hindrance.  Next,  we  determined  the  dis- 
tance r"*"  in  HCWOj  which  predicts  the  same  value  of  I].  This  HONO^  transition 
state  separation,  r'*' , corresponds  to  a specific  temperature,  given  by 
equation  (4).  An  actual  hindrance,  11^,  was  used  to  fit  the  HONO^  data  at 
this  temperature.  We  used  this  same  t/  to  predict  the  CIONO^  rates.  (This 
procedure  must,  of  course,  be  repeated  for  each  ClONOj,  temperature,  i.e.,  r"*".) 
The  larger  ClONOj  r'*'(weaker  bond,  looser  transition  state)  and  the  repulsion 
of  the  more  bulky  chlorine  atom  generally  balanced  one  another,  and  the 
hindrances  for  reactions  (1)  and  (2)  at  similar  temperatures  are  close.  High 
pressure  data  is  required  for  a more  accurate  determination  of 

Finally,  we  must  choose  an  effective,  active,  external  moment  of  inertia 
for  the  ClONOg  molecule,  a problem  not  evident  for  HONOg  because  OH,  unlike 
OCl,  could  be  approximated  as  a point  mass.  If  6 is  the  angle  between  the 
axis  of  the  active  rotation  and  the  principal  axis  , the  effective  active 
moment  is  given  by  I = I^  cos*®  + IgSin*®.  The  effective  adiabatic  moment 
(2D)  will  then  equal  the  total  ClONOj  moment  (3D)  divided  by  the  effective 
active  moment  (ID).  Two  logical  choices  for  the  active  axis  are:  (1)  An 
axis  through  the  center  of  mass  and  parallel  to  the  dissociating  bond,  or 
(2)  an  axis  parallel  to  the  active,  In-plane,  Cl-O  axis  of  rotation.  This 


second  choice  Insures  a strict  correlation  between  the  active  external 


i 

i 

I 
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molecular  rotation  and  the  active  frattment  rotations  in  the  transition  state. 

We  have  Investigated  both  possibilities.  Note  that  now  I‘*'/l  does  not  equal 
tr+  /r„)-  . 

Klgure  -1  shows  the  percent  deviation  from  third-order  low-pressure 
behavior  by  the  effective  blmolecular  rate  constant  for  recombination,  as 
a function  of  pressure.  The  results  are  Insensitive  to  the  choice  of  active 
rotor  model.  Furthermore,  the  degree  of  fall-off  Is  relatively  unchanged 
by  use  of  a different  colllslonal  efficiency.  Using  Leu,  Lin,  and  DeMore’s  value 
lk„  = 7.0  X 10"*  ‘ cm'*  s-‘  ; = l.!))"^^  at  2110",  Ak 'k^  at  50  torr  only 

Increases  by  1..  Of  course.  Inaccuracies  In  the  low-pressure  rate  constant 
will  be  reflected  In  kind  at  hlglier  pressures,  but  the  extent  of  fall-off 
will  not  be  appreciably  affected.  The  fall-off  Is  somewhat  sensitive  to  1). 

A significantly  lower  value  of  H than  that  used  Is  unlikely,  since  k Is 

CO 

already  comparable  to  the  larg<-’  HONO^  values.  Higher  1]  means  a lower  k^ 
and  a lower  pressure  commencement  of  the  fall-off.  Thus,  for  'I  = 70'’l  at 
220"  , Ak/k^,  at  50  torr  rises  by  ~ O'”!-.  Hence  the  rates  derived  from  our 
fall-off  model  can  be  viewed  as  tlie  maximum  possible  rates.  Note  in  com- 
parison to  reaction  tl),  the  degree  of  fall-off  is  smaller,  but  still 
significant  at  lower  stratospheric  pressures.  The  lower  bond  energy  of 
C lO-NO^,  means  a greater  specific  rate  constant  for  decomposition  (lesser 
for  recombination!  so  a higher  pressure  Is  required  to  insure  all  activated 
molecules  will  suffer  a stablli/.lng  collision  before  refragmenting.  Zellner-’- 
has  also  recently  done  fnlloff  calculations  for  CIONO,,  . Althougli  he  used  a 
different  value  for  k^^,,  the  degree  of  falloff  was  similar  to  the  current 
results . 

Figure  3 Includes  a similar  rough  rate  vs.  altitude  plot  for  reaction  (2! 
based  on  ovir  model  results. 
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CCWCLUSIONS 


Our  calculations  indicate  one  should  take  the  fall-off  into  account 
in  stratospheric  modeling  calculations  concerning  the  role  of  reaction  (2) 
in  trapping  active  chlorine,  which  may  destroy  the  stratospheric  layer  of 
protective  ozone.  However,  one  should  not  use  the  fall-off  of 
reaction  (1),  as  the  NAS  report  did.*^  We  believe  our  rates  are  reasonable, 
maximum  values  (minimum  fall-off,  maximum  k^) , and  thus  would  predict  the 
maximum  healing  effect  of  reaction  (2) . Certainly  future  higher  pressure 
rate  constants  for  reaction  (2)  will  test  this  point,  and  the  ability  to 
use  such  models  with  statistical  theories  to  predict  rate  constants. 

Perhaps,  future  experiments  and  modified  Gorin  model  RRKM  fits  will 
someday  permit  a more  accurate  and  systematic  method  for  the  a priori 
selection  of  hindrance  values.  The  results  here  Indicate  that  the  potential 
surface  for  HONOg  decomposition  cannot  be  fit  with  a single  parameter.  This 
has  also  been  shown  in  reference  8c,  and  future  work  may  illustrate  such 
behavior  for  a variety  of  decompositions. 
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K It'll.' RK  CART  IONS 

1.  Effective  blmoleculnr  rate  constants  for  nitric  acid  recombination 
versus  pressure  at  2St>‘’K . Experimental  points  are  from  the  data  of 
references  2 (•)  , 3 (^)  , and  -l  f ).  A wall  reaction  was  assumed 
for  the  reference  A data  The  solid  line  is  the  (3  .-18  RRKM  fit 

(parameters  in  Table  l>.  The  dashed  line  represents  termolecular 
behavior,  characteristic  of  low  pressures. 

3.  Plots  similar  to  Elgure  1 for  other  temperatures.  Data  Is  from 
reference  2 and  parameters  of  the  RRKM  calculations  are  ^jiven  In 
Table  I . 

3.  Rates  for  reactions  (1)  and  I-)  as  a function  of  the  atmospheric 

temperature  and  pres.sure  profile.  Oerlved  by  extrapolation  and  or 
interpolation  of  kjiown  experimental  data  using;  our  RRKM  calculations. 

■I  . Percentage  decline  of  the  viredlcted  RRKM  rate  constant  from  a value 
extrapolated  from  the  low-pressure  limit  for  three  calculations 
described  in  Table  1.  This  extent  of  fall-off  from  termolecular 
behavior  Is  represented  In  Elgure  1 by  the  gap  between  the  solid 
and  dashed  lines. 
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Appendix 


PARAMETERS  FOR  NITRIC  ACID  GlE'*')  AND  N(E*)  CALCULATION 


Molecule;  N(E*) 


Vibrational  frequencies;  3650,  1710,  1330(2),  880,  760, 

680,  580,  465 

Moments  of  inertia  (external):  Active  - 64.3  x 10“^° 

Inactive  - 97.3  x lO"""  (2) 


Transition  State;  G(E^) 


Vibrational  frequencies;  OH  part  - 3730 


NOg  part  - 1620,  1320,  750 


Moments  of  inertia:  Active 


OH  part  - 1.45  X 10"“°  (2) 
.92  X 10"“°  (2) 

NOj  part  - 64.3  X 10"“° 

15.45  X 10-“°(2) 
9.77  X 10-“°(2) 

Inactive,  external,  97.3  x 10"“° 


for  0%  hindrance , or 
for  60%  hindrance 

(=  external  rotation) 
for  0%  hindrance,  or 
for  60%  hindrance 

X I’*‘/I  (2) 
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Chapter  6 


A MODIFIED  CMDRIN  TRANSITION  STATE 
IN  BOND  SCISSION  REACTIONS. 

AN  APPLICATION  TO  fiTHANE  DISSOCIATION."^ 


* 
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ABSTRACT 


RRKM  calculations  are  presented  for  the  thermal  dissociation  of  ethane 
and  the  reverse  recombination  reaction  of  methyl  radicals.  The  transition 
state  employed  includes  several  hindered  rotations,  and  matches  the  exper- 
imental high-pressure  data  between  300°K  and  1400°K  using  a single  temper- 
ature dependent  hindrance  parameter.  The  experimentally  observed  rate 
constants  in  the  pressure  falloff  are  also  calculated  using  values  for  the 
average  energy  transferred  in  deactivating  collisions , In  general  accord 
with  previous  results.  The  modified  Gorin  model  is  also  consistent  with 
chemical  activation  results,  provided  the  methylene  singlet-triplet  splitting 
is  less  than  9 kcal/mole.  Comparisons  are  made  with  other  theoretical 
attempts  to  characterize  the  ethane  system. 
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I INTRODUCTION 


In  order  to  quantitatively  understand  the  mechanism  for  a complex  thermal 
chemical  process,  it  is  necessary  to  be  able  to  describe  the  elementary  thermal 
rate  constants  as  a function  of  temperature  and  pressure.  In  general,  it  has 
been  found  that  the  framework  offered  by  the  simple  transition  state  theory^  is 
sufficient.  This  generalization  is  found  wanting  only  for  simple  bond  scission 
reactions  leading  to  free  radicals  in  which  electronic  rearrangement  is  minimal. 
Thus,  in  the  system  CjHg  ^ 2CH3  to  be  discussed  herein,  a simple  fixed  trans- 
ition state  will  not  describe  the  temperature  dependence  of  the  rate  data 
adequately.  We  thus  choose  a transition  state  including  several  hindered  rot- 
ations, with  the  structure  of  the  transition  state  and  consequently  the  extent 
of  hindrance  varying  with  temperature. 

Using  the  thermod3mamic  notation  to  which  transition  state  theory  lends 
itsilf,  the  problem  may  be  simply  restated  by  pointing  out  that  the  molecular 
frequencies  necessarily  assigned  to  match  the  data  in  simple  bond  scission 
processes  require  values  of  which  in  turn  predict  large  increases  in  both 
As  (A-factor)  and  AH  (activation  energy)  over  easily  accessible  temperature 
ranges  (ca.  300-1300°K),  in  contrast  to  observations.  Quack  and  Troe®  have 
reviewed  many  simpler  systems  where  fixed  vibrational  transition  states  are  also 
Inadequate.  However,  there  is  very  little  indication  of  such  a problem  for  reactions 
characterized  by  extensive  electronic  rearrangement  in  the  incipient  radicals, 
such  as,  but-l-ene  ^ CH3  + allyl,^  ethylbenzene  ^ CH3  + benzyl,"'  and  hexa-1,4- 
diene  ^ 2 allyl.®  These  systems  showing  extensive  electronic  rearrangement  in 
the  incipient  radical  may  be  characterized  by  fixed  transition  states  with 
diminished  heat  capacity  as  a result  of  the  increasing  restriction  to  internal 
rotations.  Thus,  transition  state  theory  does  not  predict  large  temperature 
variations  in  these  cases,  in  accord  with  most  observations. 
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Simple  bond  scission  processes  occur  often,  and  we  attempt  herein  to  extend 
the  simple  approach  which  we  have  recently  utilized  in  dealing  with  the  reactions 
HO  + NO2  ^ It  • 2,®  CIO  + NO2  ClONOj,®  HOj  + NO2  ^ HOONO2,’  and  certain  ion- 
molecule  reactions  involving  the  formation  of  proton-bound  dimers  of  amines,® 

H2O,  and  H2S,®  to  the  title  reaction.  Our  future  plans  include  the  extension 
to  other  alkane  decompositions,  considering  the  importance  of  these  processes 
to  the  combustion  and  pyrolysis  of  hydrocarbons,  as  well  as  to  ion-molecule 
reactions  involving  species  of  in  erest  in  atmospheric  and  interstellar  chemistry. 

The  relevance  of  nonthermal  experiments,  such  as  chemical  activation , ‘ * 
to  the  model  presented  here  will  be  discussed  along  with  some  comments  on  current 
uncertainties  in  and  ^^(CH2L*AJ)  . 


II  BACKGROUND 

The  considerable  body  of  experimental  results  on  the  ethane  system  has 
been  reviewed  by  Waage  and  Rabinovitch**  with  more  recent  results  summarized 
by  Glanzer,  Quack  and  Troe.^®  Thus  the  related  recombination  and  decomposition 
reactions  are  good  candidates  for  theoretical  study , as  simple  proto- 
types of  Important  processes  in  hydrocarbon  combustion.  Such  attempts  provide 
a probe  of  current  assumptions  regarding  the  behavior  of  highly  excited 
polyatomic  molecules,  and  an  extrapolation  technique  for  predicting  rates  under 
different  conditions  or  for  similar  systems. 

High  pressure  data  is  available  over  a wide  temperature  range.  There  is 
considerable  spread  in  the  reported  values  for  the  high-pressure  recombination 
rate  constant  at  300-400°K,  but  the  recent  results  of  Parkes,  Paul,  and  Quinn‘® 
confirm  the  concensus  value  of  ~ 4 x 10”‘^  cm®  s“^ . The  experiments  of 
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liidleato  a sli^^htly  lower  rate  at  1UU)“K.  Kthaue  ileeomp*>si  t ton 


rates  have  been  measured  at  800-10()0‘’K  by  l,ln  and  Hack,** 
Quinn, and  can  be  Inverted  via  the  eiiul  lilirlum  constant 
ra  t «'s  . 


*”  and  by  Clark  and 
to  K i ve  reeoml>l  lui  t I on 


Several  other  pertinent  observations  have  been  nuide  on  tlie  ethane  system, 
doth  Clark  aiul  Quinn"**’  and  tilan/.er  et  al.,'^  Iiave  also  measured  rale 
constants  l'i>r  pcrtleut  eroetliane . Kxperlmental  uncertainties,  however,  preclude 
tlie  use  of  the  temperature  dependence  of  tile  isotopic  rate  ratio  as  a sensitive 
test  of  theoretical  calculations.  Ketuieily  and  Krey“’  and  Simons  and  coworkers** 
hav>'  studied  the  pressure  ilependence  of  ethane  decomposition  prepared  by  singlet 
m*'th\lene  insertion  Into  methane  (chemical  activation).  Any  model  which  fils  the 
thermal  recomlii  na  t ion  data  .should  hopt'fully  also  be  able  to  predict  these  rates  n i veil 
a ^;ood  value  for  i!)H(CH,j  ( *A , ) ) . Klnally,  pressure  dependent  rate  studies  in  the 
falloff  region  have  been  reported,*"  many  of  which  are  old  (and  scattered)  oi- 
over  a limited  ran^e  of  pressures.  Kor  a representative  test  of  our  model  over 
a wide  temperature  and  pressure  range,  we  have  considered  the  experimental 
pressure  dependence  reported  by  Gllin/.er  et  al  . *^  (l  lOO^K)  , l,ln  and  Hack*** 
(decomposition  {)ld'’K)  , and  Casas  et  al.*‘  (STO^K)  . 


Several  efforts  have  been  made  to  fit  these  results  with  IIRKM  calculations, 
as  summari.'.ed  by  Waage  and  Kablnovl  t ch ' - and  llase.**’  The  choice  of  a proper  crit  i 
configuration  is  crucial,  and  most  authors  have  chosen  a fixed  vibrational  tran- 
sition stall'  to  attempt  fitting  a wide  range  of  data.  It  is  now  clear*"’”'^  that  only 
onlv  variable  transition  state  models  fit  the  observed  temperature  dependence. 

Within  the  generally  accepted  statistical  assumption  Impost'd  by  practicalii. 
two  broad  approaches  to  the  high-pressure  rate  constant  are  possible:  the 
statistical  adiabatic  channel  theory  of  Quack  and  Troe*"*  which  counts  the  number 
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ot  uidlviilual  open  reaction  channels  whose  nuixlma  lie  at  various  positions 

alont;  a specified  reaction  coordinate;  or  alternatively,  our“  canonical  i 

approach  to  transition  state  theory  with  provisions  for  a tighter  critical 

configuration  at  higher  temperatures.  We  reason  that  the  centrifugal  iiuijtimum  closely 

describes  the  C-C  distance  in  the  transition  state  and  that  since  this  distance 

decreases  with  increasing  temperature,^-  the  transition  state  characteristics 

are  also  temperature  dependent.  We  describe  the  hindered  Gorin-^  model  which 

emerges  in  terms  of  a hindrance  parameter  which  measures  the  effective 

"tightness"  of  the  transition  state  as  a fimction  of  C-C  tiistance  and  thus, 

i 

temperature  (see  details  below). 

Quack  and  Troe‘*'have  applied  their  statistical  adiabatic  channel  (tvXC) 
iiHvlel  to  the  ethane  system,  successfully  reproducing  the  temperature  dependence 
>1  the  high-pressure  rate  constant,  although  the  absolute  numbers  are  slight  l.\ 
too  high.  They  have  not  discussed  pressure  dependence  or  chemical  activation. 

(and  they'*^)  feel  that  SAC  will  always  remain  too  difficult  to  apply  to  large  mo  1 - 
1 .•  lies  of  practical  interest.  Quack  and  Troe  have  suggested*'*^  a simpler  empi  . i < .i  1 
interpolation  technique  to  avoid  this  problem.  They  suggest  that  a one-paramt t » t 
extrapolation  from  partition  functions  of  reactant  to  those  of  products  will 
suffice.  In  the  molecules  they  consider,  the  parameter  which  they  call  ^ has 
a single  value,  but  this  is  surely  questionable  in  the  face  of  data  on  systems 

leading  to  stabilised  radicals,  such  as  mentioned  in  the  Introduction.  I'hus,  I 

these  two  ::anonical  approaches  will  probably  turn  out  to  be  very  similar  in  the 

sense  that  both  require  a system-independent  interpolation  parameter  to  maximise 

free  energy.  (A  canonical  model  never  really  addresses  itself  to  the  minimisation 

of  state  density  which  is  a mic rocanonlca 1 requirement;  thus,  the  free  energy 

chosen  might  not  corresponil  to  the  "best"  critical  configuration,  but  the 

computed  rate  constants  are  not  very  sensitive  to  minor  structural  variations 

of  the  activated  complex.)  Theoretical  predictions  in  the  falloff  region  also 

require  equal  consideration  of  the  low-pressure  rate  constant,  a factor  too  often 

given  only  cursory  attention.  Troe'-"*  and  van  den  Uergh-'’  have  applied  Tree's 

recent  low-pressure  theory, which  predicts  a collislonal  deactivation  efficiency 

(3)  from  the  average  energy  transferred,  <AK^  , to  ethane.  This  formulation  | 

will  be  followed  here. 
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III  DETAILS  OF  THE  CALCULATIONS 


The  RRKM  Model 


The  relevant  reactions  for  the  thermal  decomposition  of  ethane  are: 

k(E+) 


C^He 


2CH, 


CjHg  + M s:  CgHg  + M 


(1) 

(2) 


where  the  excited  ethane  molecules  are  those  in  the  Boltzmann  tall  above 
the  dissociation  limit  (Eq^  • According  to  the  statistical  RRKM  theory,*® 
the  decomposition  rate  is  given  by: 


k 

uni 


Q 


+^-E,/kr 

hQiQz 


P(E‘^)e 


-E'*’/kT 


dE'' 


1 + Q,+P(E+)/hQ,N(E*)F0x 


(3) 


wher  E = E — E^  is  the  maximum  energy  of  the  critical  configuration.  The 
sum  of  the  states  of  the  complex  below  energy  E'*’,  PCE"*")  , depends  on  the 
structure  chosen  for  this  transition  state.  Here,  N(E*)  is  the  density  of 
molecular  states,  Qj  is  the  partition  function  of  the  active  molecular  modes, 
Qj''‘/Qj  is  the  partition  function  ratio  of  the  inactive  modes  (I^/I,  the 
moments  of  the  adiabatic  external  rotors) , and  F is  the  Waage-Rabinovltch 
centrifugal  correction  term'*  for  conservation  of  angular  momentum  for  the 
adiabatic  rotations.  In  using  the  ethane  bond  energy  at  0°K  for  E^ , we 
effectively  assume  no  intrinsic  barrier  for  the  recombination  step  (—  1) . 

The  collision  frequency  ^ of  excited  ethane  with  the  bath  gas  M at  pressure  P 
(reaction  2)  is  calculated  from  the  Lennard-Jones  collision  cross  section  at 
temperature  T.  The  collislonal  efficiency  for  stabilization,  3,  is  calculable 
from  the  average  energy  transferred  (exponential  model)  in  reaction  (2) , 

< AE  >,  by  Troe's  formula:® 


3/(1  -/5)  = <AE>/FEkT 


(4) 
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where  the  usually  small  Fg  term  which  corrects  for  the  nonlinearity  of  the 
vibrational  state  density  near  is  given  by-* 


P _ V - 1)’.  I kT  V 

E 1=0  (s  - 1 - i)’  VEo  + / 

where  E^,  is  the  zero  point  energy  and  the  empirical  or  ~ 1 . The  recombination 
rate  constant,  k^.,  can  be  computed  from  via  the  equilibrium  constant, 

Furthermore,  the  dependence  on  E^  cancels  out. 

The  Transition  State 


For  our  model  (see  Table  I)  of  the  critical  configuration,  we  choose  the 
top  of  the  centrifugal  barrier  for  a Lennard-Jones  potential  as  the  separation 
distance  along  the  bond  axis: 

= l'*'/!  = (r*‘/ro)^  = (dOg/TlT)  ‘ (6) 

where  D = + lihv(complex)  — (molecule). 

A modified  Gorin  model® is  used  for  the  transition  state.  The  Internal 
modes  are  the  vibrations  and  rotations  of  two  methyl  radicals.  The  six 
transition  state  internal  rotations  correlate  to  four  molecular  rocking  modes, 
the  torsional  vibration,  and  the  ethane  external  rotation  about  the  C-C  bond 
axis.  This  thus  leaves  only  two  external  ethane  rotations  inactive  and 
adiabatic.  These  rotational  modes  do  not  share  in  the  random  distribution 
of  molecular  energy.  [This  treatment  conserves  the  angular  momentum  of  the 
Inactive  rotors,  but  does  not  strictly  conserve  the  total  angular  momentum. 
Note,  however,  that  the  active  external  rotor  has  the  largest  rotational 
spacing  of  the  three,  and  hence  the  lowest  average  quanta  of  angular  momentum.] 

The  four  methyl  transition  state  rotors  corresponding  to  ethane  rocking 
vibrations  are  not  actually  free  to  rotate  at  the  C-C  distances  described 
above.  The  modified  Gorin  model  hence  represents  these  degrees  of 
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Table  I 


Ik. 


PARAMETERS 


4.4 
3 .5 

(5.3)  estimated 
5.3 
3 .8 

1.31  X 10^  cm-=* 

1 .07  X 10®  cm“® 
7.75  X 10^  cra"^* 

Molecule : 

V = 2977(4),  2954,  2896,  1469(4),  1383(2), 

1190(2),  995,  822(2),  260  cm"* 

I = 11.7  X 10“^°  gm  cm* 

Inactive  I = 40  x 10“*° (2)  gm  cm* 

Transition  State:*® 

V = 3184(4),  3002(2),  1383(4),  580(2)  cm"* 

I = 11.7  X 10-*°,  29.3  X 10-*°  gm  cm* 

0*  hindrance:  29.3  x 10”*° (4)  gm  cm* 
50%  hindrance:  20.7  x 10”*° (4)  gm  cm* 


Lennard-Jones  Collision  Diameters  (A)  :** 

Ar 

I-C3H7-CO-CH3 

CH. 


freedom  as  hindered  rotations.  This  is  perhaps  a more  reasonable  picture  for 
the  "loose”  transition  states  typical  of  simple  bond  scission  reactions 

than  the  usual  very  low-frequency  vibrational  model.  This  hindrance  is 
accomplished  in  the  calculation  by  decreasing  the  effective  moments  of  inertia, 
Ijj,  of  the  two  two-dimensional  methyl  rotors; 

^ (7b) 

where  7)  is  the  hindrance  parameter  and  partition  function  for  the 

00 

active  modes  in  the  transition  state,  to  which  the  high  pressure  is 

proportional.  This  effectively  decreases  the  number  of  available  rotational 
states  by  excluding  each  rotating  radical  from  the  volume  occupied  by  the 
other,  and  thus  T1  should  rise  with  temperature,  as  r^  drops  (see  equation  6) . 
[we  note  that  in  one  sense  the  interfragment  distance  alone  does  not  entirely 
represent  the  reaction  coordinate  (r'*'^  q'*')  , since  the  four  other  modes  which 
vary  with  temperature  may  also  contribute.  However,  for  purposes  of  the  RRKM 
calculation,  the  C-C  bond  stretch  is  used  as  the  critical  mode.] 

Ethane  frequencies  and  moments  of  inertia  from  the  NSRDS  table*’  were 

used,  while  methyl  radical  values  were  taken  from  the  JANAF  tables.*®  From 

pbotoionizatlon  measurements*® of  ^f(CH3) , Eg  = 87.7  kcal/mole  was  chosen. 

From  this  information,  K was  calculated  (matching  tabulated  values)'*'*’*®  and 

eq 

Q,,  Qj  , P(E‘*')  , and  N(E*)  were  computed.  The  Lennard-Jones  cross  section  on 
which  uu  depends  was  derived  from  ethane  viscosity  data,**  and  suitably  averaged 
with  like  values  for  the  bath  gases.**  Table  I lists  some  of  the  parameters 
for  the  calculations.  Thus,  we  are  left  with  two  adjustable  parameters  at  each 
temperature  to  fit  the  wide  range  of  data.  While  < AE  >,  which  determines  the 
low-pressure  (step  2)  rate  constant,  k^ , is  expected  to  be  temperature  indepen- 
dent, T]  which  determines  the  high-pressure  (step  1)  rate  constant  k^,  is  not, 
as  mentioned  earlier.  (Thus,  if  'H  = a + bT,  we  would  have  three  constants  to 
characterize  the  system.)  In  any  event,  these  two  parameters  have  simple 
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physical  meanings,  and  any  data  fitting  attempt  is  restricted  to  a range  of 
< AE  > and  I]  values  reasonable  for  the  system. 


Chemical  Activation 

In  chemical  activation  studies,  C^Hg  is  prepared  nearly  monoenergetically 

by  insertion  of  CH2(‘Aj)  into  CH^ . The  measured  quantity,  to  be  matched 

theoretically,  is  the  ratio  of  decomposition  to  stabilization  rates,  i.e.,  the 

relative  yields  of  products  from  reactions  (1)  and  (2),  D/S.  Since  the  rate 

of  (2),  i.e.,  S, is  pressure  dependent,  the  "dimensionless"  quantity  by  which 

results  are  reported,  analogous  to  , is  u)D/S,  where  uu  is  calculated 

(roughly)  by  the  previously  described  method.  If  we  consider  deactivation  by 

a stepladder  model,  the  number  of  steps  n = (E  — Eq)/<Ae>,  D in  the  high 

pressure  limit  is  just  an  average  of  !«(£■*■)  over  n steps,  and  (uuD/S)^  = 
n 

S k(E^) . In  the  low-pressure  limit,  according  to  Robinson  and  Holbrook,^®, 
i=i 

n 

(u)D/S)g  = n k{E|)  . In  the  falloff,  the  fraction  collisionally  de- 

^ 

excited  to  the  next  step,  (u/[u)  + k(E^)J,  must  be  calculated  at  each  step,  and 

' ^ / + \ 

(juD/S  can  be  derived  from  S/(D  + -S)  fu)/(uj  + k(E^))].  We  note  that  k(Et) 
is  easily  determined  from  the  RRKM  computation,  k(E'*’)  = P(E^)/hN(E^)  ; The 
parameters  for  such  a calculation  are  U),  calculated  from  Lennard-Jones  viscosity 
diameters;  <Ae>,  determined  by  fitting  the  ethane  decomposition  data;  and  E , 
which  depends  on  the  value  chosen  for  , and  on  the  amount  of  CHjj 

vibrational  energy  from  the  photolytic  generation  not  removed  before  insertion. 
Simon's  and  Curry's  value®®  gives  E ^ 101.4  kcal/mole,  from  the  measurements 
of  McCulloh  and  Dibeler®®^  and  Zlttel  et  al.,®®  E*  iS  115.5  kcal/mole. 

A small  correction  term  must  be  added  to  the  transition  state  energy:®® 

E^  = AH®®®  — E_  + AE„„. . The  activated  ethane  after  formation  has  E = RT  ~ .6 
R o rot  rot 

kcal/mole  average  rotational  energy  in  the  two  adiabatic  external  rotors.  Since 
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angular  momentum  must  be  conserved  going  to  the  transition  state,  “ 

and  = RTtl  — — ^)  ~ .5  kcal/mole  becomes  available  to  the  other 

rot  rot  K 

modes  in  the  transition  state.  This  correction  for  average  angular  momentum 
is  akin  to  the  Waage-Rabinovi tch  correction'*'*  employed  in  the  thermal  case. 

The  fact  that  we  have  not  chosen  the  critical  configuration  in  such  a way 
as  to  assure  that  it  corresponds  to  a minimum  in  the  density  of  states  (NfE')) 
will  introduce  some  uncertainty  here.  However,  the  uncertainties  in  all  the 
input  factors  are  considerably  greater.  One  can  vary  the  critical  configuration 
in  such  a way  as  to  maintain  the  same  value  of  while  varying  AS^  (l.e., 

state  density)  and  within  reasonable  bounds. 


IV  RESULTS  .-IND  DISCUSSION 

High-Pressure  Rate  Constants  (Hindrance) 

Figure  1 shows  some  experimental  values  of  the  high-pressure  limit  of 
the  bimolecular  rate  constant  for  the  recombination  of  methyl  radicals  k^ 
as  a function  of  temperature.  The  points  between  800  and  lOOO^K  are  derived 
from  ethane  decomposition  using  the  equilibrium  constant.*®’**  The  results 
show  considerable  scatter  and  would  have  shown  more  had  we  included  some  older 
measurements.  Also  illustrated  are  the  predictions  of  statistical  adiabatic 
channel  theory,**  which  appear  slightly  too  high,  and  our  choice  of  a rough, 
empirical  fit  to  the  data.  Our  line  follows  the  perceived  slight  decline  of 
k^  with  temperature. 

This  fit  of  the  data  corresponds  to  a modified  Gorin  model  transition 
state,  as  outlined  previously,  with  a hindrance  varying  from  at  300“K 
to  81%  at  lAOO^K,  Furthermore,  this  hindrance  varies  linearly  with  r^ , the 
radical  separation  distance  in  the  transition  state,  given  as  a function  of 
temperature  by  equation  (6).  Hence  we  have 

T1  = 144  - ITr'*’  = 144  - 212  T“‘  ® i8) 


and  the  entire  fit  of  the  data  is  accomplished  by  two  parameters 


within  a physically  simple  model.  Interestingly  enough  , the 

two  constants  of  equation  (8)  can  be  related  through  the  constraint  that 
= 100"  at  r'*’  = - .6  which  is  just  twice  the  nonbonded  H-atom  interaction 
distance**  ; that  is,  the  flat  methyl  radicals  cannot  rotate  at  all  at  this 
distance.  This  would  mean  that  u at  any  one  temperature  is  sufficient  for 
predicting  ”,  at  any  other  temperature.  (A  value  of  ' near  100"  would 
indicate  that  the  four  degrees  of  freedom  which  we  are  treating  as  rotations 
should  be  treated  as  vibrations.) 

By  the  same  token  this  simple  van  der  Waals  model  would  indicate  that 
= 0"  for  r"*"  > 4.8  %,  where  the  methyl  radicals  fail  to  repel  each  other 
even  in  a C-H-  - - H-C  arrangement  11. 1 C-H  distance  + 2 x 1 .3  X H van  der 
Waals  radius  1.1  H-C  distance  = 4.8  %)  . Yet  at  300'^K  where  r’’  = 4.8  , 

”1  = 63",  not  0",  Thus,  as  in  the  nitric  acid  case  previously  explored."  the 
hindrance  is  larger  than  might  be  expected  and  more  slowly  decreasing  with 
r'*’ , although  the  ethane  values  are  more  reasonable  and  systematic.  These 
higher  hindrance  values  indicate  that  additional,  attractive  forces  are 
Important  in  restricting  or  coupling  the  motions  of  the  rotors.  The  differences 
between  the  behavior  of  T1  in  nitric  acid  and  ethane  suggest  that  the  nature  of 
the  potential  surfaces  at  these  large  separations  vary  between  systems  with 
significant  effects.  Clearly,  a larger  body  of  data  must  be  fit  with  this 
modified  Gorin  model  before  it  is  justified  by  the  observation  of  systematic 
behavior  and  available  as  a predictive  tool. 

If  ve  accept  the  values  of  C-C  dist.ince  which  correspond  to  the  minimum 
in  state  density,  from  reference  16,  we  may  calculate  i , which  is  the  value 
of  required  to  reproduce  k^  considering  these  svimewhat  smaller  distances. 

We  find  that  at  40C  K,  = 56"  compared  with  ” = 66"  and  at  860‘’K,  " * 73'. 

compared  with  "i  = 76".  Changes  of  this  nature  make  no  significant  change  in 
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kujji  over  practical  ranges  of  Interest  and  account  for  our  decision  to  forego 
the  complicated  minimization  calculation  which  they  would  require. 

The  values  of  are  somewhat  closer  to  the  van  der  Waals  model.  It 
is  possible  that  future  work  may  show  the  use  of  a microcanonical  para- 
meter, a function  of  r or  E,  with  the  minimum  density  of  states  criterion 
could  prove  more  comfortable  physically.  Our  T|  is  a canonical  parameter,  a 
function  of  r"*"  or  T,  and  our  model  an  easy  to  calculate  canonical  one. 

Shaw^*  has  recently  compared  the  forward  and  reverse  rate  constants  over 
a wide  temperature  range  for  the  reaction  CHj  + ::  CH4  + H.  The  methyl 

radical  heat  of  formation  can  be  derived  from  the  temperature  dependence  of 
the  equilibrium  constant  Kj.  = kj./kj^.  This  leads  to  an  ethane  bond  energy  of 
85.3  kcal/mole  (0°K)  , 2.4  kcal  below  the  generally  accepted  value.  In  our 
model,  T1  is  fl.xed  by  the  recombination  data  at  300°K  and  1400°K  and  is 
Independent  of  the  bond  energy.  Tj  (913°K)  shoi>ld  lie  between  1]  OGO^K)  and 
ri  (1400°K)  . Thus,  the  ethane  decomposition  data*®  at  913”-1000°K  provides 
an  added  test  of  the  bond  energy  via  our  RRKM  model.  Using  the  lower  bond 
energy  (Instead  of  87.7  kcal/mole  as  in  Figure  1)  in  the  ethane  equilibrium 
constant,  the  lOOG^K  methyl  recombination  high-pressure  rate  constant  obtained 
by  inverting  the  measured  ethane  decomposition  rate  constant  is  7 x 1G“**  cm^/s. 
This  can  be  accommodated  within  our  smoothly  varying  model  only  if  (1)  a combined 
SG'o  error  exists  in  the  recombination  (300®K)  and  dissociation  (1GG0°K)  rates. 
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which  is  unlikely  considering  the  extent  of  work  done  on  these  systems,  or 
(2)  the  1400°K  recombination  measurement‘s  is  incorrect  and  k rises  with 
temperature  (i.e.,  the  hindrance  increases  less  strongly  with  temperature). 

We  do  note,  however,  that  a 1 .5  kcal/mole  drop  in  the  ethane  bond  energy  only 
raises  the  1000°K  recombination  rate  constants  to  4 x 10”“  cm^/s  which  is 
accommodated  by  our  model,  given  reasonably  generous  estimates  of  experimental 
errors.  An  error  this  size  in  the  ethane  bond  energy  derived  from  the  dis- 
sociative photoionization  threshold  measurements^* might  be  accounted  for 
by  kinetic  energy  release  for  — CHj'*’  + H near  the  threshold.*®  The  value 

of  86.7  +.8  kcal/mole  for  DoCCHj-CHj)  derived  from  Benson  and  coworkers'  recent 
measurements*^  on  the  equilibrium  CH3  + HCl  CH^  + Cl  is  also  consistent  with 
the  Lin  and  Back  ethane  decomposition-methyl  radical  recombination  results. 
Finally,  the  838®K  ethane  decomposition  rate  constant  of  Clark  and  Quinn*”  lies 
closer  to  our  fit,  and  consequently  is  not  consistent  with  significantly  lower 
ethane  bond  energies.  (A  change  of  1.2  kcal/mole  doubles  the  recombination 
rate  constant.) 


Pressure  Dependence  (<  AK  >) 


We  have  used  an  RRKM  program^*  to  fit  the  pressure  dependence  of  the 
rate  constant  for  several  recent,  representative  experiments  which  spanned 
a large  pressure  range.  Here  ’1  is  fixed  by  equation  (8),  and  the  colllsional 
efficiency  which  determines  the  low-pressure  rate  constant  , is  calcu- 
lated from  < AE  > by  Troe's  formulation.**  Values  of  < AE  > consistent  with 
previous  determinations  would  thus  indicate  a match  of  theory  and  experiment. 

The  data  and  fits  in  Figure  2 and  the  parameters  in  Table  II  are  given 
for  the  1400®K  recombination  experiment  in  Argon  by  Glanzer  et  al.,**  the 
913“K  decomposition  results  of  Lin  and  Back  in  ethane,*®  and  the  370K  recom- 
bination results  of  Casas  et  al.  in  Isopropylmethyl  ketone  and  other  gases.** 
Excellent  fits,  suitable  for  extrapolation  with  care  to  other  pressures,  are 
obtained.  Van  den  Bergh**  has  done  similar  fits  for  the  data  of  Glanzer  et  al.* 
The  energy  transfer  values  are,  in  general,  in  accord  with  past  values*®  of 
.2  to  .7  kcal  mole"*  in  Ar  and  1.2  to  4.4  kcal  mole”*  in  C^Hg , but  the  anom- 
alously low  value  in  isopropylmethyl  ketone  (370“K)  indicates  possible 
experimental  problems.  The  theoretical  framework  provides  a good  diagnostic 
of  experimental  results. 


Table  H 

MODEL  PARAMETERS  FOR  ETHANE  DECOMPOSITION 
(Energies  in  kcal/mole) 


T(®K) 

r+(X) 

I^/I 

< AE  > 

logliv/s"*  ] 

1 

Ea 

300 

4.8 

9 .6 

63 

17  .12 

89  .56 

370 

4 .6 

9 .0 

66 

1 .0 

.41 

17.23 

89  .89 

450 

4.45 

8 .4 

67 

17  .37 

90  .16 

910 

3 .95 

6.6 

77 

2.3 

.38 

17,12 

90  .13 

1400 

3 .7 

5.7 

82 

.3(Ar) 

.06 

16.69 

88  .57 

Chemical  Activation 


Growcock  et  al.“  and  Frey  and  Kennedy‘“  have  performed  chemical 
activation  studies  on  ethane  decomposition.  The  results  give  values  for 
aJD/S  of  5.1  X 10**  and  7 .9  x 10**,  respectively.  To  calculate  theoretical 
values  of  ac/s  by  RRKM  theory  according  to  the  formulas  of  section  II,  we 
need  values  for  E^,  <AE>,  uo,  and  klE"*")  . The  Lennard-Jones  collision 
frequency,  oi,  was  calculated  for  the  experimental  gas  mixtures  from  the 
parameters  listed  in  Table  I.  The  experimental  values  of  u»/S  given  above 
are  based  on  these  values  of  la,  which  differ  slightly  from  those  used  by 
the  authors.  Two  different  values  were  considered  for  the  e.xcitation  energy 
E . These  values  are  minima,  since  CH^  vibrational  excitation  Just  prior 
to  Insertion  was  assumed  to  be  low.  We  note,  however,  that  Growcock  et  al.‘* 
claim  a 13  kcal  vibrational  excitation.  A low  value  of  99  kcal/mole  for 

■H- 

AHjCCHjj  ) taken  from  Simon  and  Curry's  “ interpretation  of  chemical  experiments 
gives  E*  2:  101.4  kcal/mole.  The  more  direct  physical  measurements  of  references 
30  and  33  combine  to  give  a high  value  of  E^  > 115,5  kcal/mole.  Welge*°  recently 
observed  singlet  methylene  from  the  laser  photolysis  (337  nm)  of  ketene.  This 
supports  a low  value  for  the  singlet-triplet  splitting  and  for  E . Recent 
theoretical  results'**  give  a splitting  of  11  kcal/mole  (E*  > 105.5  kcal/mole). 

The  most  critical  parameter  in  calculating  lUD/S  is  <AE>,  since  the  amount  of 
energy  transferred  per  collision  determines  the  number  of  collisions,  , necessary 
for  stabilization,  i.e.,  the  time  available  for  decomposition.  Since  the  bath 
gas  in  the  chemical  activation  experiments  was  predominantly  ethane,  we  have 
used  the  value  of  2.3  kcal/mole  for  < AE  > derived  previously  from  the  fit  of 
the  falloff  data  for  thermal  ethane  decomposition  at  913‘*K.  For  chemical 
activation,  however,  the  bath  gas  temperature  is  lower  (300‘*K)  and  the  excit- 
ation energies  of  the  activated  molecules  are  much  higher.  There  is  little 
definitive  evidence  on  how  < AE  > varies  with  either  of  these  parameters  .■**  •*** 


Secondly,  we  have  chosen  the  formalism  of  Troe*^  which  regards  < AE  > as  a 
constant,  over  that  of  Tardy  and  Rabinovitch**  which  considers  < AE  > down, 
the  average  energy  transferred  in  deactivating  collisions,  as  the  important 
parameter.  Since  the  Troe  method**  appears  to  be  independent  of  the  energy 
transfer  model  used,  < AE  > derived  from  an  exponential  model  for  the  SIS^K 
thermal  system  can  be  transferred  to  a stepladder  model  for  the  chemical 
activation  system.  Furthermore,  since  < AE  > » RT,<AE>  ^ < AE  > . One 

might  expect  that  any  variation  of  < AE  > with  excitation  level  would  produce 
larger  < Ae  > at  higher  E . This  would  Increase  theoretical  estimates  of  the 
stabilization  rates  S and  decrease  a®/s.  For  the  two  values  of  E*,  n = 6 
and  12 . 

Values  of  kiE'*')  are  computed  directly  by  our  RRKM  program,  with 
T = 300°K  and  Tl  = 63%,  as  given  by  equation  (8)  and  as  shown  in  Figure  3. 

They  are  calculated  before  application  of  the  rotational  energy  correction 
discussed  in  section  II , and  are  compared  with  the  J = 0 statistical  adiabatic 
channel  values.**  Our  model  predicts  slightly  lower  values  of  klE*")  at  the 
lower  E"*"  values  sampled  most  in  thermal  systems  (giving  lower  thermal  » 

but  higher  k(E^)  at  the  high  E*"  values  sampled  via  chemical  activation.  If 
it  is  the  modified  Gorin  model  which  improperly  perceives  the  potential 
surface  and  critical  configuration  (one  hopes  future  sophisticated  experiments 
will  tell),  then  the  decomposition  rate  D,  hence  the  theoretical  luD/S  must 
once  again  be  reduced. 

The  theory  predicts  minimum  high  pressure  values  for  uuO/S  of  1 .6  x 10* 
s“‘  and  4.5  x 10“*  s“*  for  the  low  and  high  values  of  AHj(CHj^)  , respectively. 
The  experiments  of  Growcock  et  al.,“  which  extend  further  into  the  high 
pressure  limit  than  Kennedy  and  Frey's  work,**  give  a value  of  5.1  x 10*  s"*  . 
Increasing  the  value  used  for  < AE  > will  lower  the  theoretical  numbers,  but 


1.54 


I 


even  in  the  strong  collision  limit,  the  high  value  for  uJ)/S  Is 

still  1 .8  X 10‘°  s~‘ . Our  theoretical  model  fits  the  experimental  result 
If  = 105  kcal/mole,  l.e.,  for  a methylene  slnglet-trlplet  splitting 


^ 9 kcal/mole  (and  our  values  of  AH^  for  CHj  and  CH^)  , In  rough  agreement 
with  the  recent  experimental"*®  and  theoretical"**  values.  Thus  our  model  is 
not  consistent  with  the  higher  value  for  AH^(CH^,  ) . 

This  is  more  clearly  Illustrated  in  Figure  4,  which  shows  several  attempted 
fits  to  Kennedy  and  Frey's  pressure-dependent  data.*®  The  deviation  from  the 
high-pressure  limit,  as  Illustrated  by  a rapidly  increasing  slope,  is  very 
apparent  In  this  range  for  the  high  value  of  E*.  Predicted  values  in  the 
falloff  region  are  quite  sensitive  to  parameter  variation.  For  example,  a 
change  In  the  value  used  for  < AE  > requires  an  almost  equal  change  in  E^. 
(However,  this  then  lowers  the  high-pressure  value  of  ojD/S.)  Our  fit  is 
consistent  with  the  similar  model  of  Kennedy  and  Frey,*®  which  used  < AE  > = 2.3 
kcal/mole,  E*  = 107  kcal/mole,  a larger  uj  and  a vibrational  transition  state 
(different  k(E)  ) . Although  the  Figure  4 falloff  fit  appears  to  be  a very 
sensitive  way  to  determine  E , an  examination  of  uncertainties  reveals  otherwise. 
Despite  our  transfer  of  the  < AE  > value  from  our  thermal  ethane  decomposition 
fit,  a 1-kcal  error  is  reasonable  and  gives  a like  uncertainty  to  E*.  Our 
differences  with  Kennedy  and  Frey  over  uJ  indicate  a 10%  uncertainty,  which  could 
change  e"*^  by  .5  kcal.  Uncertainty  in  the  position  of  the  high-pressure  exper- 
imental intercept  in  Figure  4 permits  an  adequate  fit  with  another  1-kcal 
variation  in  E . The  use  of  a vibrational  frequency  for  the  ethane  torsional 
motion  results  in  an  underestimation"*®  of  a£>/S  of  up  to  25%.  Thus  the  data 
and  our  fit  tentatively  support  the  lower  value  for  the  methylene  singlet-triplet 
splitting,  and  a range  of  higher  values,  but  not  the  19.5  kcal  mole”*  value  of 
reference  33.  While  this  does  not  consider  the  possibility  of  a much  larger 
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< AE  > at  high  excitation  energies,  or  possible  overestimation  of  k(E)  by  the 
Gorin  model,  a combined  error  of  ~ 4007o  is  required  to  support  the  high  value 
of  AH£(CH2  ).  Clearly,  however,  without  a conclusive  determination  of  E , 
this  particular  system  cannot  utilize  the  pressure  dependence  of  ud/s  to  test 
the  RRKM  theory,  k(E) , or  examine  behavior  of  < AE  > at  higher  excitation 
energies . 

The  authors  of  SAC  theory  did  not  attempt  to  fit  the  chemical  activation 
data.  We  note,  however,  that  k(E , J = 0)  is  an  upper  llmlt^^  to  k(E,J). 

Since  SAC  theory  predicts  lower  k(E)  values  than  the  Gorin  model  at  the  higher 
energies  sampled  by  chemical  activation,  SAC  values  of  UO/S  will  be  lower. 

Thus,  for  E*  = 105  kcal/mole,  E — Eq  =17  kcal/mole,  the  Gorin  value  for  k(E) 
corresponds  to  the  SAC  value  of  k(E for  an  energy  E ' of  22  kcal/mole  above 
the  decomposition  limit  (see  Figure  3).  Thus,  a rough  SAC  fit  of  the  chemical 
activation  results  requires  a higher  value  for  the  CHj  singlet  triplet  splitting, 
approximately  14  kcal/mole.  This  is  higher  than  recent  values ^ assuming 
little  CHj  vibrational  excitation,  but  well  within  all  the  uncertainties  dis- 
cussed above. 

V Conclusions 

I 

Hase,^®  and  Olson  and  Gardiner have  recently  attempted  to  predict  the 
ethane  data  with  a minimum  density  of  states  criteria  for  locating  the  critical 


configuration,  using  other  temperature-dependent  transition  state  models.  These  1 

models  portray  the  hindered  transition  state  rotors  as  methyl  rocking  vibrations  \ 

whose  frequencies  vary  with  temperature  (or  r"*")  or  as  free  rotors  with  potential  | 

■ I 

barriers  to  rotation  whose  barrier  height  varies  with  temperature.  The  simple  i 

functional  forms  used  by  the  authors  for  the  frequency  or  barrier  variation  j 

fall  to  fit  the  experimental  methyl  radical  rate  constant  temperature  variation.  j 

) 
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Physically,  more  complex  models  for  the  parameter  variations  would  be  required 
to  fit  the  data.  Furthermore,  since  the  barrier-to-internal-rotation  model 
represents  the  internal  rotations  as  free  rotors  at  the  high  energies  sampled 
by  chemical  activation,  this  theory  predicts  significantly  higher  k(E)  and 
D/S  values^’  than  our  hindered  Gorin  model.  To  summarize,  the  vibration, 
barrier,  and  Gorin  models  represent  the  transition  state  methyl  rocking  mode 
by  parabolic,  cosine,  and  square  well^®  potentials,  respectively,  and  it  is 
still  uncertain  which  form  and  resulting  k(£)  is  the  most  accurate  and  easily 
applied  one.^® 

A brief  comparison  of  the  hindered  Gorin  RRKM  model  with  the  one-parameter 
partition  function  Interpolation  procedure^"*  adopted  by  Quack  and  Troe  as  an 
approximation  to  their  statistical  adiabatic  channel  (SAC)  theory  is  of  interest. 
Noting  the  difference  in  the  k(E)  plot  (Figure  3) , any  temperature  dependent 
transition  states  constructed  to  reproduce  the  results  of  reference  14  must 
differ  from  the  hindered  Gorin  one.  Alternatively,  any  potential  surface 
interpolation  between  reactant  and  product  state  densities  which  reproduces 
the  Gorin  model  k(E)'s  must  differ  somewhat  from  the  interpolation  chosen  by 
Quack  and  Troe  to  approximate  the  SAC  results.  The  correct  transition  state, 
the  correct  interpolation,  and  which  formalism  is  more  convenient  and 
applicable,  remain  undecided. 

The  use  of  a hindered  Gorin  transition  state  in  this  work  successfully 
fits  the  data  for  ethane  decomposition  and  methyl  recombination,  using  values 
for  T]  and  AE  consistent  with  the  physical  model  and  previous  results.  These 
fits  provide  a convenient  method  for  extrapolating  the  rates  to  other,  unmeasured, 
temperatures  and  pressures.  In  contrast  to  most  vibrational  models,  the  tran- 
sition state  here  is  characterized  by  a single  temperature-dependent  parameter, 


hindrance.  The  large  entropies  required  of  the  hindered  rotations  would 
correspond  to  vibrational  frequencies  below  100  ciii~^ , and  treatment  of  the 
temperature  dependence  would  be  harder  to  visualize.  The  convenient  hindered 
Gorin  representation  is  physically  reasonable  and  meaningful  for  loose  tran- 
sition states,  and  provides  a compact  method  for  codifying  data.  More 
experience  with  the  model,  however,  is  necessary  before  a priori  predictions 
of  the  hindrance  and  estimates  of  rate  constants  become  practicable.  We  plan 
to  extend  this  work  to  other  alkane  systems  in  the  near  future. 

Secondly,  we  note  the  considerable  uncertainties  upon  which  these 
theoretical  endeavors  rest.  Besides  the  scatter  of  experimental  values, 
thermochemical  uncertainties,  and  approximate  collision  cross  sections, 

< AE  >,  and  its  temperature  variation  can  only  be  loosely  estimated.  Furthermore, 
it  is  not  established  that  the  Gorin-RRKM  model,  or  any  other,  gives  a true 
representation  of  the  sum  of  states  at  the  critical  configuration.  Any  con- 
clusions to  be  drawn  from  interpretation  of  the  chemical  activation  results 
remain  shrouded  in  the  previously  discussed  uncertainties. 

Further  research,  aimed  at  a more  microscopic  level  than  the  general 
systems  fit  here,  should  provide  more  accurate  parameters  for  such  calculations, 
test  the  Gorin  model  form  for  k(E) , and  determine  the  validity  of  the  statis- 
tical assumption  upon  which  RRKM  and  transition  state  theory  depend. 
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FIGURE  CAPTIONS 


1 High-Pressure  Methyl  Radical  Recombination  Rate  Constants  at  Various 
Temperatures 

[ i Experiments  of  Reference  18;  • Reference  13;  A Derived  from  the 
ethane  dissociation  experiments  of  Reference  19  and  Reference  20 , using 

the  equilibrium  constants  in  Table  I;  SAC  Theory  of  Reference  14;  

modified  Gorin  model  of  this  work.  The  ♦points  represent  the  inverted 
ethane  dissociation  data  for  a bond  energy  0°  = 86.2  kcal/mole. 

2 Pressure  Dependence  of  Rate  Constants 
Theory  using  parameters  of  Tables  I and  II . 

• Experiment:  (A)  Ratio  of  bimolecular  methyl  radical  recombination 

rate  constant  to  the  high-pressure  rate  constant 
at  370°K  vs.  isopropylmethyl  ketone  pressure  from 
Reference  21 . 

(B)  Ethane  decomposition  unimolecular  rate  constant 
at  913°K  vs.  ethane  pressure  from  Reference  19. 

(C)  Bimolecular  methyl  radical  recombination  rate 
constant  at  1400°K  vs.  Ar  pressure  from  Reference  13. 

3 Ethane  Decomposition  Rate  k(E)  as  a Function  of  Excess  Molecular  Energy  £~£o . 

X is  the  300°K  transition  state  of  the  modified  Gorin  model  with  no 
correction  for  angular  momentum;  • is  the  J = 0 statistical  adiabatic 
channel  result  of  Reference  14 . 


4  Pressure  Dependence  of  the  Relative  Yields  for  Reactions  (1)  (decomposition) 
and  (2) (stabilization)  for  Chemical  Activation.  Points  are  the  data  of 
Reference  10.  Ordinate  scale  is  equivalent  to  D/S  times  an  experimental 
intercept  value.  The  theoretical  curves  are  from  our  stepladder  deactivation, 
modified  Gorin  transition  state  model,  using  the  following  parameters: 

(A)  E*  = 101.4  kcal/mole  (low  CH^  singlet-triplet  splitting), 

< AE  > = 2 .3  kcal/mole; 

(B)  E*  = 115.5  kcal/mole  (high  splitting),  <^>  = 2.3  kcal/mole; 
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We  have  devised  experimental  approaches  to  the  study  of  both  the  thermochemistry 
and  thermal  kinetics  of  free  radical  reations.  In  the  course  of  this  contract 
we  have  established  thermochemical  values  for  acetonyl  radicals,  allyl  radicals, 
and  benzyl  radicals  using  our  newly  developed  techniques  of  modulated  molecular 
beam  sampling,  very  low-pressure  pyrolysis.  The  same  techniques  have  also  yield* 
rate  constants  for  allyl  combination  and  for  allyl  and  benzyl  with  HI. 
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We  have  begun  a study  of  infrared  multiphoton  chemistry  from  the  prospective  of 
gas  kineticists  interested  in  testing  theories  of  reactivity  using  non-thermal 
energisation  techniques.  In  a complimentary  study,  we  have  studied  the  thermal 
reactions  relevent  to  our  study  of  the  laser-induced  decomposition  of  ethylvinyl 
ether . 

We  have  developed  and  applied  models  for  understanding  radical-radical  inter- 
actions, as  well  as  ion-molecule  processes.  Our  use  of  a modified  Gorin  model 
has  enabled  us  to  understand  the  pressure  and  temperature  dependence  of  these 
processes . 
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